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ABSTRACT 

Context. Understanding the relation of the attenuation with the gas is fundamental for interpreting the appearance of galaxies. This 
can now be done down to a local scale in the local Universe, thanks to the high spatial resolution achievable from the FUV to the FIR. 
More importantly, this relation is also crucial for predicting the emission of galaxies. This is essential in semi-analytic models that 
link dark matter from large cosmological simulations, to the baryonic matter which we can observe directly. 

Aims. The aim of the present paper is to provide new and more detailed relations at the kpc scale between the gas surface density and 
the face-on optical depth directly calibrated on galaxies, in order to compute the attenuation not only for semi-analytic models but 
also observationally as new and upcoming radio observatories are able to trace gas ever farther in the Universe. 
Methods. We have selected a sample of 4 nearby resolved galaxies and a sample of 27 unresolved galaxies from the Herschel 
Reference Survey and the Very Nearby Galaxies Survey, for which we have a large set of multi-wavelength data from the FUV to 
the FIR including metallicity gradients for resolved galaxies, along with radio HI and CO observations. For each pixel in resolved 
galaxies and for each galaxy in the unresolved sample, we compute the face-on optical depth from the attenuation determined with 
the CIGALE Spectral Energy Distribution fitting code and an assumed geometry. We determine the gas surface density from HI and 
CO observations with a metallicity-dependent X C o factor. 

Results. We provide new, simple to use, relations to determine the face-on optical depth from the gas surface density, taking the 
metallicity into account, which proves to be crucial for a proper estimate. The method used to determine the gas surface density or the 
face-on optical depth has little impact on the relations except for galaxies that have an inclination over 50°. Finally, we provide detailed 
instructions on how to compute the attenuation practically from the gas surface density taking into account possible information on 
the metallicity. 

Conclusions. Examination of the influence of these new relations on simulated FUV and IR luminosity functions shows a clear impact 
compared to older oft-used relations, which in turn could affect the conclusions drawn from studies based on large scale cosmological 
simulations. 

Key words, galaxies: ISM, ISM: dust, extinction, Infrared: galaxies, Radio lines: galaxies, Ultraviolet: galaxies 



1. Introduction 

Understanding galaxy formation and evolution is one of the cen- 
tral questions in modern astrophysics. Major efforts in terms 
of observations and theory have been made over the last few 
decades in order to gain insight into the physical processes driv- 
ing the transformation of the protogalaxies seen at high redshift 
into the diverse range of galaxies we see in the local Universe. 
Large scale numerical simulations have become increasingly 
important tools over the last 10 years with remarkable break- 
throughs. With recent advances in computing power, cosmolog- 
ical simulations, which were initially limited to massive struc- 
tures, now have a mass resolution sufficient to resolve objects 



down to dwarf galaxies in large volumes (e.g. Guo et al. 2011) 
yielding unprecedented insight into the formation and evolution 
of galaxies across a large dynamical range of types and masses. 

What determines the formation and evolution of galaxies is 
the transformation of baryonic matter. Yet the computational 
cost of an extensive physical model of the processes affecting 
baryonic matter (gas heating and cooling, star formation, feed- 
back, active nuclei, metal enrichment, etc.) at the sub-grid level 
remains particularly high. Even large hydrodynamical simula- 
tions, such as the Overwhelmingly Large Simulations project 
(Schaye et al. 2010), that model in detail the physics of the 
baryons are either limited to significantly smaller volumes or 
require a dazzling amount of computing power: 10 5 cores and 
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65.5 x 10 9 particles in the case of the MassiveBlack simulation 
(Di Matteo et al. 2012) for instance. Large scale cosmological 
simulations such as the Millennium simulation (Springel et al. 
2005), which are limited to dark matter, are then coerced to em- 
ploy coarser, yet physically sound, "recipes". Still, the baryonic 
matter is key to constrain the formation and evolution of galax- 
ies as it is directly traced through high redshift surveys. The 
connection between the dark matter traced by N-body simula- 
tions and the baryonic one traced by observations can be made 
through SAMs (semi-analytic models) such as the ones devel- 
oped by Croton et al. (2006) or De Lucia et al. (2006) for in- 
stance, among many others. Reassuringly, comparisons between 
different SAMs show that they are globally consistent but some 
systematic uncertainties can nevertheless be found (Fontanot 
et al. 2009a). To constrain models, simulated observations can be 
computed from the combination of numerical simulations with 
SAMs. These artificial catalogues can then be compared with 
observations (e.g. Henriques et al. 2012). 

The attenuation due to the dust affects dramatically 
the observed or predicted fluxes, colours, and Spectral 
Energy Distributions (SEDs) of galaxies which are reddened. 
Computing the attenuation in galaxies is an important chal- 
lenge observationally, and one of the most delicate parts in 
SAMs (Fontanot et al. 2009b) along with the computation of the 
amount of cold gas (Obreschkow et al. 2009). As the radiation 
from different stellar populations propagates through clouds of 
gas and dust, it finds itself partly absorbed and reprocessed to- 
wards longer wavelengths. Because the energetic radiation from 
the youngest stellar populations is the most affected by the pres- 
ence of dust, the most widely used star formation tracers are 
strongly impacted, either positively (the infrared) or negatively 
(the ultraviolet or Ha). This heavily affects our ability to trace 
star formation and by extension our understanding of galaxy for- 
mation and evolution across cosmic times. 

In galaxies, the observed attenuation is actually the conse- 
quence of an intricate interplay between 1) the relative weights 
of the different stellar populations across the SED, 2) the relative 
geometry of the dust and the stars, and 3) the line-of-sight of the 
observation. Radiation transfer codes can be used to model in 
detail the emerging light of dusty galaxies (Baes & Dejonghe 
2001; Tuffs et al. 2004; Bianchi 2008; Popescu et al. 2011). 
Unfortunately their high computational cost makes them unsuit- 
able for cosmological simulations. In SAMs, different methods 
to compute the attenuation yield different predicted counts in 
ultraviolet (UV) Luminosity Functions (LFs) by a factor ~ 2 
(Finlator et al. 2006; Fontanot et al. 2009b). This is worrisome 
as the UV is one of the most widely used star formation trac- 
ers (Kennicutt 1998; Kennicutt & Evans 2012) and is observable 
by ground-based surveys for galaxies with z > 1. One of the 
methods followed to compute the attenuation in SAMs revolves 
around variations of the following: 1) link the face-on optical 
depth to the gas surface density which is a known quantity, and 
2) assuming a given extinction curve, geometry, and line-of- 
sight, retrieve the V-band attenuation from the face-on optical 
depth (e.g. Guiderdoni & Rocca-Volmerange 1987; Devriendt 
et al. 1999; De Lucia & Blaizot 2007). Even though galaxies 
have a complex geometry, such methods are generally calibrated 
on objects in the Milky Way. In our galaxy, the reddening can 
be linked with the gas column density (Burstein & Heiles 1978, 
1982; Schlegel et al. 1998). However, studies carried out by Xu 
et al. (1997) and Boissier et al. (2004, 2007) on different samples 
of galaxies actually found only weak relations between the face- 
on optical depth and the gas surface density. This sheds doubt on 
the reliability of such methods, and rises the issue of the poten- 



tial impact on our understanding of galaxy formation and evo- 
lution. The main difficulties in performing studies on galaxies 
essentially revolved until recently around the coarse resolution 
available to trace both the gas and the attenuation in galaxies. 
The increasing availability of high resolution maps of the gas 
through HI and CO observations, and multi-wavelengths obser- 
vations from the FUV to the FIR now allow us to gain a much 
needed insight into the relation between the gas and the atten- 
uation in galaxies. Beyond the direct impact of such relations 
on SAMs, with the inception of ALMA and other millimetre in- 
struments, we can trace the gas content of galaxies across the 
Universe both indirectly through the emission of the dust and di- 
rectly through molecular lines such as CO. The gas can then be 
used to correct the SED for the attenuation. Therefore it is now 
timely to reexamine at last the relation between the attenuation 
and the gas surface density in galaxies. 

The overarching aim of this paper is to determine whether 
there is a physical relation between the face-on optical depth 
and the gas surface density in galaxies at kpc scales and if so, to 
provide new, well calibrated, but still simple and easy to use rela- 
tions to compute the attenuation for observations and for models. 
To do so we carry out a study on 2 samples made of nearby re- 
solved and unresolved galaxies defined in Sect. 2. The sensitive 
and delicate issue of the determination of the attenuation and of 
the gas mass surface density is detailed in Sect. 3. We provide 
a first overview of the relation between the attenuation and the 
difference gas phases in Sect. 4. We then derive new relations 
between the face-on optical depth and the gas mass surface den- 
sity for resolved galaxies (Sect. 5). In Sect. 7 we describe how 
to apply these new relations to obtain a reliable estimate of the 
attenuation. We then examine the case and unresolved galaxies 
(Sect. 6). We compare the new relations we have obtained m to 
standard relations from the literature in Sect. 8. The impact of 
the new relations on LFs is discussed in Sect. 9, and finally we 
conclude in Sect. 10. 



2. Sample and data 

2.1. Sample selection 

2.1.1. Multi-wavelength requirements 

We want to determine the relation between the attenuation and 
the gas in nearby galaxies. To reach this goal we need to deter- 
mine these two quantities on a resolved basis in galaxies as well 
as for unresolved galaxies. 

The attenuation which is due to absorption and scattering of 
the radiation out and into the line-of-sight can be determined 
in several ways. A widespread method to correct for the atten- 
uation in nearby galaxies is to compute the ionised, gas-phase 
attenuation from the B aimer decrement between the Ha and H/3 
lines (e.g. Lequeux et al. 1979; Kennicutt 1992, and many oth- 
ers). This is best performed using spectra, and would ideally re- 
quire integral field unit spectroscopy of the disk to map the at- 
tenuation reliably when considering resolved galaxies. Narrow- 
band observations of galaxies can be adversely affected by un- 
derlying absorption features, especially for the H/3 line. In any 
case this method requires stellar populations no older than typ- 
ically 10 Myr. Another method is to combine extinguished and 
unextinguished star formation tracers such as UV or Ha with 
24 yum for instance. However, such methods often rely on a sta- 
tistical analysis over a sample of galaxies (Kennicutt et al. 2009; 
Calzetti et al. 2010) or a sample of star forming regions (Calzetti 
et al. 2005, 2007; Kennicutt et al. 2007; Relano et al. 2007; Li 
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et al. 2010) rather than being specifically tailored to each indi- 
vidual region. With the availability of large multi-wavelength 
datasets from the FUV to the FIR at sufficient resolution to re- 
solve nearby galaxies, it is now possible to determine the attenu- 
ation accurately from the continuum emission by modelling the 
SED not only of entire galaxies but also on a pixel-by-pixel ba- 
sis, naturally taking into account the local physical parameters 
(e.g. Boquien et al. 2012). Ly et al. (2012) showed that at least 
in the case of unresolved galaxies, the use of SED fitting to deter- 
mine the attenuation is as robust as using the Balmer decrement 
from spectra. We therefore determine the attenuation from a full 
SED modelling. 

The gas is overwhelmingly made of hydrogen either in 
atomic or molecular form. The atomic gas is easily detectable 
through its 21 cm spin-flip transition. The detection of molec- 
ular gas however is more difficult due to the intrinsic nature of 
H2, which accounts for the bulk the molecular gas mass. The 
most widespread method is to observe a rotational transition of 
the CO molecule which is a good tracer of molecular gas even 
though it is dependent on the temperature or the density of the 
gas (Maloney & Black 1988), and its metallicity (Wilson 1995; 
Israel 1997, 2000; Barone et al. 2000; Israel et al. 2003; Strong 
et al. 2004; Israel 2005). 

Following the multi-wavelength requirements we have ex- 
posed, in order to examine the relation between the gas and the 
attenuation at local and global scales in galaxies, we need 1) high 
resolution HI and CO data to trace the atomic and molecular gas, 
and 2) high resolution multi-wavelength data from the FUV to 
the FIR to determine the attenuation. For unresolved galaxies the 
wavelengths requirements are the same. 

2.1 .2. Targets and data 

We carry out the current study in the context of the Herschel 
SPIRE guaranteed time consortium. We rely on 2 surveys dedi- 
cated to the observation of galaxies in the nearby Universe: the 
Herschel Reference Survey (HRS, Boselli et al. 2010) and the 
Very Nearby Galaxies Survey 1 (VNGS). We use these surveys 
as a baseline to guide our selection of a resolved sample (here- 
after sample A), and an unresolved sample (hereafter sample B). 

For sample A we select all late-type galaxies from these 2 
surveys that are not dominated by a strong active nucleus, and 
that have high resolution data available in all bands at a sufficient 
depth: UV, optical, NIR, and FIR to model the SED accurately 
and determine the local physical parameters (including the atten- 
uation), and also HI and CO to trace the atomic and molecular 
gas. To decide whether a galaxy is sufficiently resolved or not, 
which depends on its size, distance, and inclination, we rely on 
visual inspection of the data at 350 fim. We drop the 500 /im 
band because its resolution is coarser and it does not add value 
to compute the attenuation. If structures are visible, the galaxy 
is deemed resolved for our purpose. We carry out the analysis at 
a resolution of 30", only slightly coarser than that of the 350 fim 
band. The shallow depth of GALEX UV data for some galaxies 
causes further reduction in the number of elements in the sam- 
ple. The final set of photometric bands is the following: FUV, 
NUV (GALEX), u', g', r', i', z' (SDSS), J, H, Ks (2MASS), 
70 /urn (Spitzer/MTPS), 250 fim, and 350 pm (Herschel/SPIRE), 
in addition to HI and CO maps. Despite the tremendous develop- 
ment of millimetre and radio imaging, high quality data cover- 
ing the entire body of nearby galaxies at a resolution better than 

1 http : //herschel . esac . esa . int/Docs/GTKP/KPGT_ 
accepted . html#KPGT_cwilsoQl_l 



30" remain rare. This is the most stringent criterion limiting the 
number of galaxies in the sample. We rely on two large surveys 
for HI data using the Very Large Array (VLA): VIVA (Chung 
et al. 2009) and THINGS (Walter et al. 2008). For CO, we rely 
on the HERACLES survey (Bigiel et al. 2008) carried out at the 
Institut de Radioastronomie Millimetrique (IRAM), which cov- 
ers a sample of nearby galaxies in the CO(2-l) transition at 6" to 
map the molecular gas. In conclusion, sample A is constituted of 
4 galaxies: NGC 2403, NGC 4254, NGC 4321, and NGC 5194 
(Table 1). 

For sample B, we have applied the same basic criteria on the 
original HRS sample, substituting the requirement on the reso- 
lution for a similar requirement on the inclination bja > 0.33 
in order to limit radiation transfer effects. We use the following 
combination of photometric bands which maximises the size of 
the sample: FUV, NUV, U, g\ V, r' , i' , J, H, Ks, 60 fim (IRAS), 
250 yum, and 350 /urn, in addition to HI and CO measurements. 
The requirement on the availability of both HI and CO is very 
stringent, along with the low-inclination requirement, both elim- 
inating a large number of candidates. Sample B is made of 27 
galaxies (Table 2). 

Finally we summarise the origin of the data used in this study 
in Table 3. 



2.2. Data processing 

The data processing requirements are different for both samples. 
We must ensure that sample A allows us to perform a pixel-by- 
pixel analysis across the entire SED. Conversely, for sample B 
we can use directly the fluxes that have been published in the 
literature so no extra processing is required. 

For galaxies in sample A, the processing of the FUV to FIR 
bands, including HI and CO maps is performed similarly to 
Boquien et al. (2012). In a nutshell: background and foreground 
objects are masked; UV and optical data are corrected for the 
foreground Galactic attenuation using the attenuation value from 
Schlegel et al. (1998) combined with a Cardelli et al. (1989) 
extinction curve including the O'Donnell (1994) update; using 
the convolution kernels published by Aniano et al. (2011), im- 
ages are convolved to a Gaussian with a 30" Full Width Half 
Maximum (from 0.5 kpc for NGC 2403, up to 2.5 kpc for 
NGC 4254 and NGC 4321) corresponding to the lowest reso- 
lution of the entire dataset; they are projected on the same refer- 
ence grid with a pixel size of 10"; the background is subtracted; 
finally pixels under a 3-<x threshold in any band, are discarded. 
The standard deviation of each image is computed taking into 
account both the variation of the background and the pixel-to- 
pixel noise. The entire procedure is described in greater detail in 
the aforementioned article. 



3. Attenuation and gas surface mass density 

3.1. Determination of the attenuation 

As mentioned in Sect. 2.1.1, to determine the attenuation, we 
model the SED for each pixel in each galaxy (sample A), or of 
unresolved galaxies (sample B). We present here how this mod- 
elling is performed, and what is the expected accuracy. The im- 
pact of using alternative methods such as the relation between 
the IR-to-UV ratio and the attenuation is discussed in appendix 
A. 
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Table 1. Sample A. The physical parameters of the galaxies are provided by Boselli et al. (2010); Ciesla et al. (2012) (see references 
therein), or obtained from the NASA Extragalactic Database. 



NGC# 


RA 


Dec 


Type 


K 


a 


b/a 


Inclination 


Distance 




(J2000) 


(J2000) 




(Vega mag) 


(kpc) 




(degrees) 


(Mpc) 


2403 


07 36 51.40 


+65 36 09.2 


SAB(s)cd; HII 


6.19 


10.7 


0.50 


60 


3.37 


4254 


12 18 49.63 


+ 14 24 59.4 


Sa(s)c 


6.93 


30.4 


0.91 


24 


17.0 


4321 


12 22 54.90 


+ 15 49 20.6 


SAB(s)bc; LINER; HII 


6.59 


45.1 


0.89 


27 


17.0 


5194 


13 59 52.71 


+47 11 42.6 


SA(s)bc pec; HIISy2.5 


5.50 


11.3 


0.83 


34 


6.92 



Table 2. Sample B. The physical parameters of the galaxies are provided by Boselli et al. (2010) (see references therein), or obtained 
from the NASA Extragalactic Database. 



HRS# 


NGC# 


RA 
(12000) 


Dec 
(J2000) 


Type 


K 

(Vega mag) 


a 

(kpc) 


b/a 


Inclination 

(degrees) 


Distance 
(Mpc) 


015 


3338 


10 42 07.54 


+ 13 44 49.2 


SA(s)c 


8.13 


29.8 


0.62 


52 


17.4 


077 


4030 


12 00 23.64 


-01 06 00.0 


SA(s)bc;HII 


7.33 


23.6 


0.72 


44 


19.5 


081 


4045 


12 02 42.26 


+01 58 36.4 


SAB(r)a;HII 


8.75 


14.8 


0.74 


42 


17.0 


084 


4067 


12 04 11.55 


+ 1051 15.8 


SA(s)b: 


9.90 


5.9 


0.74 


42 


17.0 


089 


4178 


12 12 46.45 


+ 10 51 57.5 


SB(rs)dm;HII 


9.58 


26.5 


0.35 


70 


17.0 


096 


4212 


12 15 39.36 


+ 13 54 05.4 


SAc:;HII 


8.38 


17.8 


0.56 


56 


17.0 


100 


4237 


12 17 11.42 


+ 15 19 26.3 


SAB(rs)bc;HII 


10.03 


10.0 


0.58 


55 


17.0 


102 


4254 


12 18 49.63 


+ 14 24 59.4 


SA(s)c 


6.93 


30.4 


0.91 


24 


17.0 


114 


4303 


12 21 54.90 


+04 28 25.1 


SAB(rs)bc;HII;Sy2 


6.84 


32.6 


0.81 


36 


17.0 


122 


4321 


12 22 54.90 


+ 15 49 20.6 


SAB(s)bc;LINER;HII 


6.59 


45.1 


0.89 


27 


17.0 


142 


4393 


12 25 25.50 


+16 28 12.0 


Sa? pec;HII 


9.49 


12.9 


0.50 


60 


17.0 


154 




12 26 47.23 


+08 53 04.6 


SAB(s)cd 


10.71 


23.1 


1.00 





23.0 


157 


4420 


12 26 58.48 


+02 29 39.7 


SB(r)bc: 


9.66 


10.0 


0.42 


65 


17.0 


171 


4451 


12 28 40.55 


+09 15 32.2 


Sbc: 


9.99 


13.2 


0.49 


61 


23.0 


196 


4519 


12 33 30.25 


+08 39 17.1 


SB(rs)d 


9.56 


17.8 


0.72 


44 


17.0 


203 


4532 


12 34 19.33 


+06 28 03.7 


IBm;HII 


9.48 


12.9 


0.38 


68 


17.0 


204 


4535 


12 34 20.31 


+08 11 51.9 


SAB(s)c;HII 


7.38 


41.2 


0.89 


27 


17.0 


205 


4536 


12 34 27.13 


+02 11 16.4 


SAB(rs)bc;HII;Sbrst 


7.52 


35.8 


0.45 


63 


17.0 


206 




12 34 39.42 


+07 09 36.0 


SBm pec;BCD 


11.01 


9.9 


0.70 


46 


17.0 


208 


4548 


12 35 26.43 


+ 14 29 46.8 


SBb(rs);LINER;Sy 


7.12 


29.7 


0.83 


34 


17.0 


221 


4580 


12 37 48.40 


+05 22 06.4 


SAB(rs)a pec 


8.77 


10.7 


0.74 


42 


17.0 


246 


4651 


12 43 42.63 


+ 16 23 36.2 


SA(rs)c;LINER 


8.03 


19.3 


0.71 


45 


17.0 


247 


4654 


12 43 56.58 


+ 13 07 36.0 


SAB(rs)cd;HII 


7.74 


24.6 


0.52 


59 


17.0 


256 


4691 


12 48 13.63 


-03 19 57.8 


(R)SB(s)0/a pec;HII 


8.54 


12.1 


0.81 


36 


14.8 


287 


4904 


13 00 58.67 


-00 01 38.8 


SB(s)cd;Sbrst 


9.50 


11.9 


0.71 


45 


17.0 


293 


5147 


13 26 19.71 


+02 06 02.7 


SB(s)dm 


9.73 


8.1 


0.81 


36 


14.5 


307 


5364 


13 56 12.00 


+05 00 52.1 


SA(rs)bc pec;HII 


7.80 


32.5 


0.65 


49 


16.5 



Table 3. Origin of the selected data. 



Galaxy/sample 


UV 


Optical 


NIR 


FIR 


HI 


CO 


NGC 2403 


GALEX/GI3-50 


SDSS 


2MASS 


SAG2/VNGS 


VLA/THINGS 


IRAM/HERACLES 


NGC 4254 


GALEX/GI2-17 






B08+SAG2/HRS+HeViCS 


VLA/VIVA 




NGC 4321 


GALEX/NGA 












NGC 5194 


GALEX/GI3-50 






SAG2/VNGS 


VLA/THINGS 




Sample B 


GALEX/GUViCS 


SDSS+SAG2/HRS 


2MASS 


IRAS+SAG2/HRS 


Arecibo 


KittPeak 12m 



References. 2MASS: Skrutskie et al. (2006); Arecibo: Boselli et al., in prep.; B08: Bendo et al. (2008); GALEX: Martin et al. (2005); GUViCS: 
Boselli et al. (201 1); Cortese et al. (2012); HeViCS: Davies et al. (2012); IRAM/HERACLES: Leroy et al. (2009); IRAS: Neugebauer et al. (1984); 
Sanders et al. (2003); Kitt Peak 12m: Boselli et al., in prep.; NGA: Gil de Paz et al. (2007); SAG2/HRS: Boselli et al. (2010); Ciesla et al. (2012); 
SAG2/VNGS+HeViCS Spitzer/MWS: Bendo et al. (2012b) SDSS: Abazajian et al. (2009); VLA/THINGS: Walter et al. (2008); VLA/VIVA: 
Chung et al. (2009). 

Notes. The resolution of each dataset is the following. GALEX: 6"; SDSS: 1.4"; 2MASS: 2.5"; Spitzer/MIPS: 18"; Herschel/SPIRE: 18-25"; 
VLA:6-25";IRAM: 6". 



3.1.1. CIGALE 

We use the SED modelling code CIGALE (Burgarella et al. 
2005; Noll et al. 2009) which handles both the emission of the 
stars in the UV-optical domain, and the absorption and reemis- 
sion of this radiation by the dust at longer wavelength through 



an energy balance requirement. Several recent works have made 
use of it to determine the physical parameters of low and high- 
redshift galaxies (Buat et al. 2011a,b; Burgarella et al. 2011; 
Giovannoli et al. 2011) as well as subregions within galaxies 
(Boquien et al. 2012). CIGALE has been amply described in the 
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aforementioned articles. We succinctly remind its main charac- 
teristics here. 

3.1.2. Physical parameters 

CIGALE handles a large number of input parameters that can be 
fine-tuned to achieve different goals. This means that we have 
to select the parameters space carefully to constrain the attenua- 
tion precisely. We use the set of parameters defined in Boquien 
et al. (2012) as a baseline. We briefly describe the choice of these 
parameters hereafter. 

The stellar component which dominates the UV, optical, and 
NIR domains is made of 2 populations: the first one modelling 
the old stellar population, and the second one modelling the 
most recent burst of star formation. Both populations follow a 
Kroupa (2001) initial mass function with a standard metallicity 
of Z = 0.02. We take the latter parameter as fixed because our 
samples are made of objects that have typically a solar neigh- 
bourhood metallicity and there is no object with an extreme 
metallicity that would induce a strong effect on the blue end of 
the stellar spectrum. The relative masses of these 2 populations 
can vary freely to take into account the strength of the last star 
formation episode compared to the older stellar population. The 
emission of the stellar populations in the UV and in the optical 
is attenuated following a bump-less starburst-like attenuation 
law (Calzetti et al. 1994, 2000) of varying slope, parameterised 
by a (A/Aq) 6 multiplying factor, with A the wavelength, A® the 
normalisation wavelength, and 6 the slope-modifying param- 
eter. Allowing for a differential attenuation between the older 
and younger stellar populations, the energy is absorbed by the 
dust following the previously defined attenuation curve and is 
re-emitted in the IR, conserving the total energy. The dust emis- 
sion is modelled using the star-forming galaxies templates from 
Dale & Helou (2002), parameterised by the IR power-law slope 
a. We should mention that because the computation of the IR 
luminosity is largely independent from the chosen library, this 
choice of the template library has no impact on our results. In 
Tab. 4, we present the list and the values of the parameters used. 

Using this space of parameters, we have generated a total of 
1455300 models with CIGALE. We present in Fig. 1 a compari- 
son between observed and modelled flux ratios constraining star 
formation, the attenuation, and the properties of the IR emission. 
We see that these models and the observations overlap, showing 
that this set of CIGALE models is adequate to reproduce the ob- 
servations. In the case of NGC 5 194, we notice that there is a tail 
towards high r' /NUV and 10/FUV flux ratios. These regions 
correspond to the northern end of the galaxy that is contami- 
nated by NGC 5195, a peculiar SB0 galaxy (Sandage & Bedke 
1994), with which NGC 5194 is interacting. 

Even if the models can reproduce the observations, it is not 
known how accurately the intrinsic physical parameters can be 
retrieved from such modelling with this combination of bands. 
Ensuring that the attenuation is reliably determined is especially 
important to understand its relation with the gas surface density. 
To do so, we follow the procedure outlined in Giovannoli et al. 
(201 1) and Boquien et al. (2012) for instance. In a first step, we 
fit the SED of each pixel in each galaxy in sample A. Using the 
best fits, we then create a catalogue of artificial SEDs and their 
associated parameters which are perfectly known. Subsequently 
to model uncertainties on the observed fluxes, we add a random 
flux to each band in each artificial SED, following a Gaussian 
distribution whose width is determined from the standard devi- 
ation of the original data. This allows us to obtain an artificial 
catalogue of simulated observations of perfectly known regions 



in galaxies that are at the same time representative of the origi- 
nal sample. Finally, we compare the parameters of the artificial 
galaxies determined through the analysis of the probability dis- 
tribution function by CIGALE, with the intrinsic, true parame- 
ters. The analysis of the results shows that the determination of 
the FUV attenuation is more reliable than the V-band attenuation 
(Fig. 2). We note however a slight underestimate of the attenu- 
ation for the highest values (Apuy >2). These biases remain no 
larger than the typical uncertainty which is computed from the 
probability function. As such, they should not affect the results 
of the paper. We see that the relative scatter of the attenuation 
in the FUV band is typically 30% to 50% smaller than for the V 
band, while the Spearman correlation coefficients are higher. For 
this reason, and because FUV more closely traces star formation 
than optical bands, we will derive relations between the gas sur- 
face density and the FUV attenuation rather than the V-band 
attenuation as is customary. We provide a systematic method to 
extend the FUV attenuation to all wavelengths in Sect. 7. 

3.2. Determination of the atomic and molecular gas mass 
surface density 

3.2.1. Sample A 

The atomic and molecular gas components can be traced indi- 
rectly through the emission of the dust, or in a more standard way 
through HI and CO observations. We adopt the latter method as a 
baseline in this paper as it allows us to ensure that the determina- 
tion of the gas surface density is entirely independent from the 
determination of the attenuation, which would not be the case 
with the former method. The impact of using the dust to trace 
the gas is discussed in appendix B. 

As long as HI remains optically thin, the conversion from HI 
emission to surface mass density is straightforward and reliable. 
For sample A, we convert the maps from brightness tempera- 
ture to gas mass surface density using the following relation (see 
appendix C for the derivation): 

8.85 

i 

with Y,hi the atomic gas surface density in M Q pc~ 2 , Aa and AS 
the major and minor axes of the beam in arcsec, S , the HI flux 
density in Jy per beam in channel i, and AV the channel width 
in m s _1 . We have to note however that self-absorption is not 
always negligible. A recent study by Braun et al. (2009) on the 
inclined galaxy M3 1 found that the HI mass had to be increased 
globally by 30% to correct for self-absorption, and even more 
locally. The lower inclination of our galaxies and the relatively 
coarser resolution should limit the problem of HI self-absorption 
and its spatial variation. 

The rotational transitions of CO are in theory good tracers of 
the molecular gas because its critical density is low enough to be 
excited collisionally in the bulk of the molecular gas. The rela- 
tion between CO emission and the amount of molecular gas has 
been extensively studied in the literature but remains subject of 
significant uncertainties that depend on the metallicity, the tem- 
perature, the density, etc. Ideally, the combination of several CO 
transitions permits to constrain the physical properties of the CO 
emitting regions to trace the molecular gas as accurately as pos- 
sible. Unfortunately such a large dataset is not available for the 
current sample. While CO(1-0) maps are often the most com- 
monly available, for sample A we use the CO(2-l) maps from 
the HERACLES survey (Leroy et al. 2009) because they provide 



5 



Boquien et al.: Towards understanding the relation between the gas and the attenuation in galaxies at kpc scales 
Table 4. List of CIGALE parameters. 



Parameter 



Unit 



Range 



Metallicity 

E-folding time of the old population 

Age of the old population 

E-folding time of the young population 

Age of the young population 

Mass fraction of the young to old stellar population 

Initial mass function 

Slope modifying parameter 6 

V-band attenuation of the young population 

Reduction of the attenuation for the old population 

IR power-law slope a (Dale & Helou 2002) 



Z 

Gyr 
Gyr 
Myr 
Myr 



mag 
mag 



0.02 
12345 67 
13 

5 10 20 50 100 200 
5 10 20 50 100 200 
0.0001 0.001 0.01 0.1 0.999 
Kroupa (2001) 
-0.4 -0.3 -0.2 -0.10.0 
0.05 0.15 0.30 0.45 0.60 0.75 0.90 1.05 1.20 1.35 1.50 
0.25 0.50 0.75 
1.0 1.5 2.0 2.5 3.0 3.5 4.0 
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Fig. 1. Comparison between observed and modelled flux ratios. The shade of the grey hexagons represents the logarithm of the 
density of models, with a darker shade indicating a higher number of models in the bin. The observed fluxes are indicated by circles, 
the colour depending on the galaxy. The median error bars including systematic uncertainties are shown in the upper-left side of 
the plots. The left plot shows the r' /NUV versus NUV/FUV flux ratios and the right plot shows the 10/FUV versus r' /NUV flux 
ratios. 
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Fig. 2. Attenuation derived from the analysis of the probability distribution function performed with CIGALE versus the true value 
for artificial regions in galaxies, in order to test the accuracy of the modelling. In each quadrant we indicate the number of data 
points N, and the Spearman correlation coefficient p at the top and the scatter of the relative difference in the attenuation at the 
bottom. The solid black line indicates the 1-to-l relation. The error bar on the left indicates the median uncertainty as computed 
from the analysis of the probability distribution function. For consistency with the analysis carried out later in this paper, we have 
only selected datapoints with a good best fit (x 1 < 2). The left panel describes the FUV band attenuation whereas the right one 
describes the V band attenuation. 
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simultaneously a good resolution and a good sensitivity. We dis- 
cuss in appendix D the impact of using different CO transitions. 
To convert the CO emission to molecular gas surface density, in 
a first approach, we will use a standard Xco conversion factor 
valid for galaxies that have a metallicity similar to that of the 
Milky Way (Strong et al. 1988): 

X co = 2.3 x 10 20 cm" 2 (K km s -1 )~\ (2) 
the molecular gas surface density coming directly: 

Z Hl = 2m p x X co x I C0 , (3) 

with m p the mass of a proton, and 7 C0 , the CO line inten- 
sity. Even though there are more recent estimates, it is well 
within the typical range of Xco found in the local group (1-4 x 

10 cm " 2 (Kkms 1 ) ). As we will see later the choice of this 
value has been made for consistency with Boselli et al. (2002). 
As mentioned before, one of the main parameters affecting the 
Xco factor is the metallicity, with the conversion factor being 
higher in less metallic environments owing to complex factors 
linked to differences in dust shielding, lower CO abundance etc. 
The impact of such a dependence on the metallicity is discussed 
in Sect. 5.3. 

The Xco factor is generally defined for the CO(l-O) line. As 
shown by Braine & Combes (1992), the typical CO(2-l)/CO(l- 
0) ratio for star forming galaxies does not show large variations 
from one galaxy to another: 0.89 ± 0.06. We therefore divide 
CO(2-l) maps by this factor before computing the molecular gas 
surface density with Eq. 3. We should note that this ratio could 
be lower in the interarm regions (in the range 0.4-0.6) compared 
to the arms (in the range 0.8-1.0) as can be seen in NGC 5194 
for instance (Koda et al. 2012). 

Finally, to account for the presence of Helium, we multiply 
the gas masses by a factor 1.38. This value is derived assuming 
X = 0.7154 and Y = 0.2703 (Asplund et al. 2009). All quoted 
gas masses will implicitly be inclusive of Helium. 

3.2.2. Sample B 

For unresolved galaxies, the computation of the mean gas sur- 
face density is more complex as it depends on the large scale 
distribution of the gas and the intrinsic characteristics of the ob- 
servations (beam size, pointing strategy, etc.). The atomic gas is 
generally distributed on a large scale length and can extend well 
beyond the optical disk, fueling star forming regions seen in the 
UV with GALEX (Thilker et al. 2005, 2007). Conversely on a 
broad scale the molecular gas is mostly concentrated towards 
the centre of the galaxy with an exponentially declining profile. 
High-resolution observations reveal however much more com- 
plex profiles with important variations in the azimuthal distri- 
bution (Regan et al. 2001). Similarly to Xu et al. (1997), we 
consider different distributions to compute the mean gas surface 
density. When considering the atomic hydrogen only, we assume 
the gas has an exponential distribution and we average over a ra- 
dius rui corresponding to 3 scale-lengths, rni — 3 x 0.61 x r2s 
(Bigiel & Blitz 2012), with r25, the optical radius of the galaxy 
at 25 mag arcsec~ 2 . When considering the molecular hydrogen, 
which is traced through CO(1-0) for all galaxies in sample B, 
we also assume an exponential distribution and we similarly av- 
erage over 3 scale-lengths, r Hl = 3 x 0.2 x (Leroy et al. 2008; 
Lisenfeld et al. 2011). When considering both we simply sum 
the atomic and molecular gas surface densities. We caution that 



these are only coarse estimates of the actual distribution of the 
gas and variations from galaxy to galaxy are to be expected. 

Similarly to sample A, all masses are corrected to take into 
account the presence of Helium. 

3.3. Other gas components 

We have considered here only atomic gas traced through the 
21 cm line and molecular gas traced through CO. However, that 
may not account for the total mass of gas in the galaxy as it does 
not take into account the "dark gas" and the ionised gas. 

3.3.1. Dark gas 

Besides the molecular gas traced by CO, there can also be so- 
called "dark gas" which is not traced by CO and therefore may 
be missed. This gas can be found in large quantities in rel- 
atively low metallicity galaxies such as the Large Magellanic 
Cloud for instance (Galliano et al. 2011), and also in a hand- 
ful of high redshift objects (Maiolino et al. 2009), perhaps also 
due to the lower metallicity of these objects compared to lo- 
cal spiral galaxies. However our sample is made of relatively 
metal-rich spiral galaxies for which CO should be a more reli- 
able tracer of the total molecular content. Still, some dark gas 
can also be found in the Milky Way (Grenier et al. 2005; Langer 
et al. 2010; Pineda et al. 2010; Velusamy et al. 2010; Planck 
Collaboration et al. 2011). However, theoretical models suggest 
that for a given mean attenuation, the fraction of dark gas re- 
mains constant (Wolfire et al. 2010). Considering this, we have 
decided not to explicitly take into account this component. If 
necessary, it can be taken into account afterwards, rescaling the 
derived relations. 

One way to naturally take into account the presence of this 
dark gas is to compute the gas surface density indirectly from 
the emission of the dust. We discuss this in detail in appendix B. 

3.3.2. Ionised gas 

The last component we discuss here is the ionised gas. It can 
be found both in the disk in regions that have recently seen 
star formation and in the warm medium surrounding galaxies 
(e.g. Ferriere 2001; Haffner et al. 2009). Even though some dust 
has been detected in the warm medium (Lagache et al. 1999, 
2000), its contribution to the optical depth is difficult to ascer- 
tain. Consequently, we choose here to consider only the gas 
ionised by young, massive stars in the disk. This may induce 
an increase in the uncertainties on the determination of the fi- 
nal attenuation. However, regions with recent star formation are 
also regions where the atomic and molecular gas have the high- 
est surface density, which makes their relative importance un- 
certain. Studies show that the emission of the warm gas is more 
closely associated to star formation than the emission of the cold 
gas (Bendo et al. 2010; Boquien et al. 201 1; Bendo et al. 2012a). 
In nearby galaxies, the warm to cold dust mass ratio has been 
shown to vary from ~ 100 to a few 1000s (e.g. Willmer et al. 
2009; Kramer et al. 2010), suggesting that only a small fraction 
of the dust actually lies in ionised gas. To ensure this, we have 
computed the approximate mass of ionised gas in each pixel. We 
have first computed the ionising photon luminosity following the 
case B recombination (Osterbrock 1989): 

Q(H°) = ^7rR 3 n 2 e a B [photons s _1 ] , (4) 
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with R the radius of a spherical HII region, n e the mean electron 
density and ag the case B recombination coefficient for the H/? 
line. Doing so we have assumed a constant density of ionised hy- 
drogen within the radius R and a fully neutral hydrogen beyond. 
It easily comes that: 



Q(H°) 



hi- 



6 A/ 

Expressing Q (h } as a function of the Hy8 luminosity: 



(5) 



UHp) 



V jnR 3 a B hv 



(6) 



with h the Planck constant and v the frequency of the Hfi line. 
Accounting for He the total gas mass associated with the HII 
region can then be easily calculated: 



M H ii = 1.38m# 



\%nR 3 L(H0) 



[Mo], 



(7) 



with niH the mass of an hydrogen atom. 

We compute the H/3 luminosity from the SFR. First we con- 
vert the SFR to an Ha luminosity. Then we consider a temper- 
ature of 10000 K, which yields L(Ha) /L (H/3) = 2.87. If we 
assumed that each 10" pixel contains typically an HII region 
with a radius of 50 pc, we find that (Lhii) = 0.3 ± 0.4 M pc~ 2 , 
with a maximum of Y.hii - 3.5 M Q pc~ 2 . While these estimates 
remain uncertain because of the number and sizes of HII regions 
within a single pixel, it shows that the total mass of the gas in HII 
regions is very small compared to the atomic and molecular gas 
components we have considered so far. Thus, we have decided 
not to account for the ionised gas mass to determine the relation 
between the attenuation and the optical depth. 
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4. A first overview of the relation between the gas 
and the attenuation in resolved galaxies 

In Fig. 3 we plot for resolved galaxies the pixel-by-pixel re- 
lation between the FUV band attenuation derived in Sect. 3.1 
and the gas surface density derived in Sect. 3.2. We note that 
there is little correlation between the attenuation and the atomic 
gas in general Qp\< 0.18), except for NGC 2403 (p = 0.72). 
Conversely, we observe mild correlations with the molecular gas 
for all galaxies (p > 0.47) except for NGC 2403 (p = 0.27). This 
difference between the atomic and molecular gas can be useful 
to understand where the attenuation takes place in galaxies as 
we will see hereafter. A good relation can also been seen when 
considering the total (HI+H2) gas surface density for all galaxies 
(p>0.47). 

Interestingly, we see that the relation between the attenu- 
ation and the gas column density follows 2 regimes. At high 
column density (typically beyond 10-20 M pc -2 ), there is a 
clear increase of the attenuation with the gas column density. 
Conversely, under this threshold, the attenuation seems to be 
rather independent from the gas column density and converges to 
a value circa 1.5 mag. This is best seen in the molecular gas, but 
can also be noted when the total gas column density is consid- 
ered. We will examine the physical origin of this dual behaviour 
in Sect. 5.2. 

At low molecular gas mass surface density, we see an in- 
crease of the scatter in Afuv- This is at least in part due to larger 



Fig. 3. Attenuation in the FUV band versus the atomic (top), 
molecular (middle), and total (bottom) gas mass surface density 
in units of M Q pc -2 for individual regions in galaxies (sample 
A). Each data point corresponds to a single pixel which has been 
detected at a 3-cr level from the FUV to the FIR. No selection 
has been applied on CO and HI data. Data points with a poor fit 
(X 2 ^ 2) have also been removed. At the bottom of each panel we 
display the number of selected pixels N and the Spearman corre- 
lation coefficient p. Finally in the bottom right corner the median 
uncertainty on Afuv computed by CIGALE is indicated. 

relative uncertainties on the determination of Apuy. Indeed, 
these regions tend to be on the outskirts of star-forming regions, 
with lower fluxes and therefore larger relative errors than in more 
active regions. This ultimately propagates to the determination 
of the FUV attenuation. 

The reason why we see a correlation between CO and dust 
attenuation in some galaxies but not others may be related to 
changes in the dust heating among galaxies. The particular case 
of NGC 2403 is interesting in that it shows an opposite be- 
haviour compared to the other elements of the sample. Bendo 
et al. (2010) found that in NGC 2403 the CO emission does 
not correlate well with tracers of cold dust when examined as 
small spatial scales, which may partly explain the absence of a 
good correlation between the molecular gas and the attenuation. 
Along this line, a recent study of Bendo et al. (2012a) has shown 
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that NGC 2403 has a thermal component that is not heated by 
the star forming regions. A similar analysis on NGC 4254 and 
NGC 4321 (Bendo et al. 2013, in prep.) does not find such a 
dust component. The recent results of Mentuch Cooper et al. 
(2012) on NGC 5194 also suggest that the dust is mainly heated 
by star forming regions. Hence, in NGC 4254, NGC 4321, and 
NGC 5194, the molecular gas and most of the dust attenuation is 
associated with the FIR emission near regions of star formation. 
In NGC 2403 however the dust attenuation associated with FIR 
emission is not in the same place as the molecular clouds, which 
causes the relation to break down. Though, since the FIR emis- 
sion in NGC 2403 is at least partly from dust absorbing light 
in the diffuse interstellar medium (ISM), since the atomic gas is 
also in the diffuse ISM and not in star forming regions, the dust 
attenuation appears correlated with atomic gas. While studying 
further the relation between the star formation, and the attenu- 
ation is of utmost interest to improve our understanding of star 
formation across the Universe, this would go well beyond the 
scope of the present paper. We refer to the aforementioned pa- 
pers for further details. We nevertheless briefly discuss the re- 
lation between the molecular gas fraction and the attenuation in 
appendix E. 

5. Relation between the gas and the face-on optical 
depth in resolved galaxies 

As we have shown earlier, in external galaxies the attenuation 
and the gas surface density are related in a non-trivial way. 
Indeed, whether a photon is affected or not by dust depends both 
on the relative geometry of the stars and the dust, but also on the 
line of sight of the observer. Unfortunately, no analytic model 
can fully grasp the impact of complex geometries of galaxies. 
They are bound to be coarse approximations. This may be worri- 
some as for a given quantity of dust some geometries will cause 
a more efficient attenuation than others, and edge-on galaxies 
will be more attenuated than face-on galaxies. If the attenuation 
is dependent on the inclination of the galaxy in a complex way, 
it is much easier to correct for it if we work on face-on optical 
depths. The overarching aim of this paper is to provide relations 
to compute the attenuation rather than estimates of the face-on 
optical depth per se, which are a means and not an end. Therefore 
we adopt a 2-step approach. First, we assume a given geometry 
to compute the relation between the face-on optical depth and 
the gas mass surface density (Sect. 5 and 6). In a second step, 
these relations can be applied on any object, assuming the same 
geometry, to retrieve the attenuation, as explained in Sect. 7. The 
use of the same geometry largely cancels out its initial impact on 
the determination of the relation between the attenuation and the 
face-on optical depth (appendix F). 

In this section we present the stars-dust geometry we adopt 
(Sect. 5.1) in order to compute the relation between the attenu- 
ation and the face-on optical depth (Sect. 5.2), and we examine 
the impact of the metallicity in Sect. 5.3. 

5.1. Geometry 

5.1 .1 . Choice of the geometry 

As described by Disney et al. (1989), several geometries be- 
tween the dust and the stars can be reasonably considered in 
galaxies, such as a simple dust screen, a sandwich model (a thin 
layer of dust embedded in a thick layer of stars), a slab (an in- 
finite plane-parallel geometry in which stars and dust are per- 
fectly mixed over the same scale), or even a clumpy ISM. While 



the screen model can be adapted for studies in the Milky Way, it 
is probably too simplistic to reproduce the complex relative dis- 
tributions of stars and dust in external galaxies. The slab model 
is more realistic as gas is intimately mixed with stars in star- 
forming regions. In the sandwich model, the relative heights of 
the dust and the stars can be freely adjusted to take into account 
their different scales (e.g. Boselli et al. 2003). Another common 
model is that of a clumpy ISM (e.g. Witt & Gordon 1996, 2000). 
The clumps of dust can be modelled following a Poisson distri- 
bution, the source being partially attenuated by the presence of a 
variable number of clumps that all have the same size and optical 
depth (Natta & Panagia 1984). 

More complex models, such as the one presented by Chariot 
& Fall (2000) for instance, beyond defining a geometry also take 
into account the differential attenuation between stellar popula- 
tions of different ages, as younger stars are still embedded in 
their dusty birth clouds whereas older ones have broken out of 
these clouds. Such models are by nature very dependent on the 
star formation history (SFH) of galaxies. However as we con- 
centrate here exclusively on the FUV band where the youngest 
stellar populations dominate, such models do not present a sig- 
nificant advantage over simpler ones. To ensure this, we have 
carried out a supplementary run of CIGALE without differen- 
tial attenuation. It yields FUV attenuations very close to when a 
differential attenuation is considered showing the FUV band is 
indeed dominated by one stellar population. 

As stated earlier, the adopted geometry is a means to correct 
for the inclination and determine intrinsic relations between the 
optical gas and the gas surface density. As long as the geometry 
is fully consistent end-to-end, its choice is not critical. To de- 
termine the relation between the face-on optical depth and the 
gas mass surface density, we have decided to adopt the widely 
used slab geometry. We examine the impact of using different 
geometries in Sect. 5.1.2 and appendix F. 

For a slab geometry, if we assume an isotropic scattering 
the attenuation and the face-on optical depth are then related 
through the following equation adapted from Natta & Panagia 
(1984): 



-2.5 log 



' 1 - exp (— VI - to) at a /cos i) 
VI - o)aTa/cos i 



(8) 



with A,i the attenuation and the face-on optical depth at wave- 
length A, i the inclination angle of the galaxy, and u>x the albedo 
at wavelength A. The isotropic scattering is introduced through 
the term yl - u>a- This is only a simple approximation. To take 
scattering fully into account radiative transfer models would be 
necessary (Baes & Dejonghe 2001; Pierini et al. 2004; Tuffs 
et al. 2004; Inoue 2005; Rocha et al. 2008), which is much be- 
yond the scope of the present article. Nevertheless, for the re- 
mainder of the paper we take cjfuv = 0.38, following the values 
published by Draine (2003b) for a Milky Way extinction curve 
withfl v = 3.1. 

To compute the face-on optical depth from the attenuation, 
we have inverted Eq. 8: 



cos i 



Ta 



VT 



^(-lO^x^^'^+lO - 4 ^], (9) 



with W the Lambert W function. The derivation of this relation 
is provided in appendix G. 
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5.1 .2. Impact of the inclination correction on the derived 
attenuation 

As we will see later, the relations we determine to compute the 
face-on optical depth from the gas surface density depend sig- 
nificantly on the geometry. The reason is that a given attenua- 
tion corresponds to very different optical depths depending on 
the geometry. This is not necessarily a problem as long as the 
geometry is kept consistent. The application of these relations 
would yield identical attenuations if all galaxies were face-on 
and if there was no scatter around the relations. As our data de- 
viate from these ideal conditions, to evaluate the actual impact 
of the geometry we present in Fig. 4 the ratio between the at- 
tenuation obtained assuming a given geometry and that obtained 
assuming a slab geometry, for the gas surface density spanned 
by sample A. We see that while using relations other than for a 
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Fig. 4. Ratio of the FUV attenuation obtained from a relation for 
a given geometry to the attenuation for a slab geometry versus 
the gas surface density. We assume an inclination of 45°, corre- 
sponding to the mean of datapoints in sample A, which is neces- 
sary to compare adequately with relations for clumpy geometries 
as the number of clumps can vary according to the inclination. 



slab geometry induces an offset, generally towards a lower at- 
tenuation, it is limited to 15% at most, in case we consider an 
extreme, screen geometry. This shows that our original assump- 
tion that the geometry would only have a limited impact on the 
results is indeed validated. In Fig. 5, we show how the relative 
attenuations vary as a function of the inclination angle and the 
gas surface density, assuming slab and screen geometries. We 
see that for inclinations no larger than 50°, we get attenuations 
within 20% from either a screen or a slab geometry. However 
the attenuations become increasingly discrepant for higher incli- 
nations. This translates both the uncertainties due to the scatter 
around the relations, and probably the complexity of the geome- 
try within galaxies compared to these simple models. We there- 
fore warn of the greater uncertainties when deriving the atten- 
uation from the face-on gas column density for highly inclined 
galaxies. 



Fig. 5. Ratio of the FUV attenuation assuming screen and slab 
geometries. The inclination angle varies from 0° (blue) to 90° 
(red). 

using Eq. 9. To compute the gas column density, in this first step 
we do not take into account the impact of the metallicity which 
will be explored in detail in Sect. 5.3. We present in Fig. 6 the 
relation between the face-on optical depth and the gas surface 
density. Combining data points from all regions of the galaxies 
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Fig. 6. Face-on optical depth in the FUV band versus the total 
gas mass surface density in units of M pc~ 2 for a slab geometry, 
and assuming a metallicity-independent X co factor. The blue 
curves represent the best fit for the full sample whereas the black 
ones represent the best fit for individual galaxies for which the 
equation is given at the top of each panel. The uncertainties on 
Tfuy are indicated by the error bar at the bottom right of each 
panel. 

in sample A, we find the following relation between the face-on 
optical depth and the gas surface density: 



5.2. Gas surface density and face-on optical depth 

Following the choice of a slab geometry in the previous section, 
we have computed the face-on optical depth in the FUV band 



t fuv = (2.226 + 0.040) + (0.058 ± 0.001) x E H , (10) 

for 2.5 < Y.H < 221.0 M pc 2 . The presence of a non-zero con- 
stant is surprising at first sight. A first explanation could be that 
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there is a selection effect against FUV-bright, low-attenuation 
regions with a low gas surface mass density, as we require emis- 
sion at a 3-cr level in all bands. However, the inspection of se- 
lected regions shows that except for the most external regions, 
the entire extent of galaxies has been selected, within the limit 
of the observed field-of-view. We conclude that this is unlikely 
to be the reason for the observed levelling at low gas mass sur- 
face density. A more likely explanation would be that this is ac- 
tually an effect due to the mixing of different regions because of 
the level of the resolution. Indeed, at a resolution of 0.5 kpc to 
2.5 kpc (30" beamwidth) there is a blending between different 
regions. This affects differently the measure of the gas mass sur- 
face density and the attenuation. The computation of the result- 
ing gas mass surface density is straightforward, being simply the 
average of all regions within the resolution element. However, 
the attenuation is strongly weighted towards the brightest re- 
gions in the FIR and in the FUV band. As a consequence, if 
we consider that a given resolution element actually contains 2 
regions that have respective gas mass surface densities S[ and S2 
covering the same area and with si » S2, the measured gas sur- 
face density over the entire resolution element will be ~ si/2. If 
we assume that these regions have attenuations ci\ and «2, if the 
first region is much brighter than the second one, the measured 
attenuation in the resolution element will the ~ a\. We there- 
fore see that such blending will differently affect the gas mass 
surface density and the attenuation. In other words, the latter is 
luminosity weighted whereas the former is governed by a simple 
geometrical mean over the surface of a pixel. This is even more 
evident if we consider that denser regions tend to form stars more 
actively, and then will drive the attenuation within the resolution 
element. The consequence is that if there is blending between 
different regions within a galaxy there will be a plateau at low 
gas mass surface density, which is what we see here. This means 
that the relations we derive in this paper should not be applied at 
a resolution more refined than that considered here. Conversely, 
this suggests that standard relations that do not take into account 
such blending should not be applied to unresolved galaxies or 
sufficiently large regions. 

The second aspect that is to be noted is that there is a 
large scatter, both within galaxies and from one galaxy to an- 
other. The internal scatter can largely be explained by the un- 
certainty on the face-on optical depth. We have to note that this 
uncertainty is computed from the uncertainty on the FUV at- 
tenuation provided by CIGALE. When comparing the relation 
provided in Eq. 10 to relations derived for individual galaxies 
(Fig. 6) there are clear differences. The slope is relatively sim- 
ilar for NGC 4254, NGC 4321, and NGC5194, but it is about 
twice as steep for NGC 2403. Conversely, NGC 2403 levels 
at a much lower face-on optical depth than average, whereas 
NGC 4254 and NGC 4321 level at a higher face-on optical 
depth. NGC 5194 is close to the value for the entire sample. This 
is probably due to intrinsic differences between galaxies in the 
sample. Another possibility would be that it is due to the physi- 
cal resolution. However, after having carried out a similar study 
at a resolution of 64" we found that such a coarsening had no 
effect on the relative levels. Incidentally this suggests that be- 
yond a certain resolution, the level of the face-on optical depth 
at low column density is left mostly unchanged. Indeed if the 
mixing within a resolution element is already representative of 
the content of the galaxy, increasing the mixing will not change 
the properties within a resolution element. 

Finally, we see that for NGC 4321, there is a visible satura- 
tion in the optical depth at high column density. This undoubt- 
edly affects the fit. Performing a fit only when < 60 M Q pc~ 2 



yields a steeper slope for NGC 2403 than when all regions are 
considered but only change marginally the parameters for other 
galaxies. 

5.3. Impact of the metallicity 

A change in the metallicity can have dramatic effects on the mea- 
surement of the molecular gas reservoir of galaxies. The reason 
is that the molecular content is generally traced through quanti- 
ties that are dependent on the metallicity. For instance the Xqo 
factor mentioned in Sect. 3.2 is expected to vary with the metal- 
licity. Ignoring this effect naturally increases uncertainties on the 
gas mass derived from the CO emission. At lower metallicities 
CO lines trace an increasingly small fraction of the total molec- 
ular gas, directly increasing the value of the Xqo factor, barring 
any other effect. 

Several studies have been dedicated to quantifying the rela- 
tion between the metallicity and the Xqo factor. Wilson (1995) in 
a seminal study found a sub-linear relation between the Xqo fac- 
tor and the oxygen abundance. Arimoto et al. (1996) and Boselli 
et al. (2002) found nearly linear slope (-1.01 and -1.00, respec- 
tively). Recently, Schruba et al. (2012) found especially steep 
slopes between -2.8 and -2.0. Conversely, Blitz et al. (2007) or 
Bolatto et al. (2008) did not find any relation with the metallicity. 

To investigate the influence of the metallicity on the rela- 
tion between the gas surface density and the face-on optical 
depth, we examine how this relation changes when we consider a 
metallicity-dependent Xqo factor. To do so we adopt the relation 
derived by Boselli et al. (2002), which is intermediate between 
the steep slopes found by Schruba et al. (2012), and the lack of 
relation found by Blitz et al. (2007) or Bolatto et al. (2008): 

\ogX co IX COsal = 8.92 - 1.01 x [12 + logO///] , (11) 

with Xco gan the standard Milky Way factor. To be consistent with 
the metallicity scale used by Boselli et al. (2002), we adopt the 
metallicity gradients measured by Zaritsky et al. (1994) rather 
than the newer gradients published by Moustakas et al. (2010). 
This prevents the introduction of biases due to differences be- 
tween metallicity estimators. In Fig. 7, we plot the relation be- 
tween the total gas surface density and the face-on optical depth 
when using the relation given in Eq. 1 1 . Computing the relation 
between the gas surface density and the face-on optical depth, 
we find: 



t fuv = (2.060 + 0.059) + (0.137 ± 0.005) x S H . (12) 

We note that there is a clear steepening of the relations with re- 
spect to those shown in Fig. 6. This is expected because Eq. 11 
naturally yields a lower (respectively higher) Xqo factor for 
12 + log O/H > 8.83 (resp. 12 + log O/H < 8.83). 

Beyond affecting the estimate of the molecular gas surface 
density, the metallicity also has a direct impact on the gas-to- 
dust mass ratio and therefore on the amount of dust attenuating 
the interstellar radiation. The simple relations we have derived 
do not fully capture the impact of the variation of the metallicity, 
in effect averaging over all the line-of-sights, objects, and metal- 
licities. To take into account these variations, different methods 
can be used. A first, very naive, method would be to correct the 
face-on optical depth computed as in Eq. 1 2 for the difference 
of the gas-to-dust mass ratios between sample A and that of the 

observed object: = r FUV x ^Z!^^ with ^uv 
the metallicity corrected face-on optical depth, Tfuv (£h) the 



11 



Boquien et al.: Towards understanding the relation between the gas and the attenuation in galaxies at kpc scales 





14 




12 




10 




8 








6 




4 




2 




14 




12 




10 




8 








6 




4 




2 









10 



NGC 2403 / 
N=1158 p = 0.72 /// 


NGC 4254 : / 
N=514p = 0.44. It 










NGC4321 t 
N=676 p = 0.49 II 


NGC 5194 / 
N=1427 p = 0.55 1 1 











9.18 



36 
30 
24 
1 

1.12 C 

1.06 c 

1.00 

1.94 

1.88 

1.82 



S H [M a pc- 2 ] 



Y, H [M 8 pc- 2 ] 



Fig. 7. Face-on optical depth in the FUV band versus the to- 
tal gas mass surface density in units of M Q pc~ 2 . The colour of 
each point indicates the oxygen abundance determined from the 
gradients published by Zaritsky et al. (1994). The black line rep- 
resents the best fit for the entire sample (Eq. 12) when taking 
into account a metallicity-dependent Xqo factor (Eq. 1 1 ). The 
coloured lines represent the best fit for the entire sample (Eq. 13) 
for different values of 12 + log O/H, when fitting simultaneously 
the gas column density and the oxygen abundance. At the top of 
each panel the number of datapoints and the Spearman correla- 
tion coefficient is indicated. 

face-on optical depth obtained from Eq. 12, and MujMdusi the 
gas-to-dust mass ratio for sample A (subscript sample A), and 
for the targeted object (subscript object). Indeed, as mentioned 
before, a variation of the gas-to-dust mass ratio will directly 
affect the quantity of dust affecting the emission. This ratio is 
evidently strongly dependent on the metallicity. Therefore if is 
important that the metallicity of the sample the relation was de- 
rived on is homogeneous and well-known. With a metallicity 
spread over more than 0.5 dex in sample A, the former crite- 
rion is certainly not fulfilled as there are significant variations of 
the gas-to-dust mass ratio from one region to another, and from 
one galaxy to another. A better way to take the metallicity into 
account is to derive relations between the face-on optical depth 
and the gas surface density, parameterising them explicitly on 
the oxygen abundance. With this method, we find: 



tfuv = [(1-926 ± 0.034) + (0.051 + 0.002) x I H ] 

x j q(0.947±0.024)([ 12+log O/H] -9.00) 



(13) 



The normalisation of the oxygen abundance, here 9.00, is arbi- 
trary. We have chosen this value because it is close to the typical 
values yielded by the oxygen abundance estimator of Zaritsky 
et al. (1994) for nearby spiral galaxies. This relation is necessar- 
ily dependent on the way the metallicity has been determined as 
different methods give different estimates. One can use the rela- 
tions provided by Kewley & Ellison (2008) to convert the metal- 
licity computed from one method to the Zaritsky et al. (1994) 
method. Note that because of the metallicity and the gas sur- 
face density are not completely independent quantities in galax- 
ies (regions with a higher gas surface density tend to be more 
metal rich than regions with a lower gas column density), great 
care must be taken to interpret the coefficients of the best fit re- 
lated to the gas surface density and the metallicity. Examining 



the variance-covariance matrix of the best fit parameters, we find 
that the correlation coefficient between the coefficients related to 
these quantities is r — -0.63. 

In Fig. 8 we show the face-on optical depth computed with 
these relations versus the face-on optical depth derived from the 
attenuation provided by CIGALE and an assumed slab geome- 
try. Taking into account the metallicity brings as expected sub- 




r KN , measured 

Fig. 8. FUV face-on optical depth computed with Eq. 12 taking 
into account the gas surface density (green dots), and Eq. 13 
taking into account the gas surface density and the metallicity 
(red dots), versus the measured face-on optical depth from the 
attenuation provided by CIGALE and an assumed geometry. The 
black line represents to 1-to-l relation. 



stantial improvements. The most notable change is due to the re- 
gions belonging to NGC 2403 which have a systematically lower 
metallicity than other galaxies in the sample. The regions be- 
longing to NGC 2403 correspond to the concentration of points 
in the lower left corner in Fig. 8. If the impact of the metallicity 
is not taken into account, Eq. 12 would systematically overes- 
timate the optical depth. Eq. 13 brings a strong improvement, 
yielding optical depths much closer to the ones obtained from 
CIGALE attenuations. 

The relative difference of the face-on optical depth goes 
from: 



\ TFUV, measured 



0.43 ± 0.77, 



with Atfuv - t fuv, fit - tfuv, meaun 
the gas surface density, to: 



e 4 when considering only 



At 



FUV 



T FUV, measured 



= 0.11+0.39, 



(15) 



when also taking into account the metallicity. This means that 
the uncertainty on the face-on optical depth is of the order of 
~ 50%. This hints at the importance of additional, undetermined, 
parameters, and it shows the inherent complexity of determining 
the attenuation, and by extension the face-on optical depth, from 
the gas column density. 

For the remaining of this article we will always consider a 
metallicity dependent Xqo factor unless specifically noted other- 
wise. 
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6. Comparison with unresolved galaxies 

Due to their more complex, large scale geometries, gas distribu- 
tion, and stellar populations, the relation between the gas surface 
density and the face-on optical depth in unresolved galaxies is 
expected to be more challenging to determine. It is important to 
verify whether the relations previously obtained provide reliable 
results for unresolved galaxies as it conditions their applicabil- 
ity. One of the first attempts to connect the attenuation and the 
gas surface density in unresolved galaxies was made by Xu et al. 
(1997). They found weak relations between the attenuation in 
the B band and the face-on optical depth. More recent studies 
have yielded similar results (Boissier et al. 2004, 2007). 

A first inspection of the relation between the FUV attenua- 
tion and the gas surface density in sample B is shown in Fig. 9. 
There is a weak anti-correlation between the attenuation and the 



In Fig. 10 we show the relation between the face-on optical 
depth and the gas surface density. In spite of the scatter and the 
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Fig. 9. Same as Fig. 3 for sample B. 

atomic gas mass surface density (p = -0.33). Conversely, there 
is a clear correlation with both the molecular (p = 0.64) and the 
total gas content (p = 0.54). This is actually due to the fact that 
some galaxies have a high molecular fraction. If attenuation is 
indeed linked to denser regions, then it is expected that galaxies 
with relatively little HI due to the high molecular fraction will 
be more attenuated. 
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Fig. 10. Relation between the FUV face-on optical depth and the 
gas surface density for unresolved galaxies (black circles with 
error bars). The blue line corresponds to the best linear fit for 
the unresolved sample. The corresponding relation is indicated 
at the top of the figure, along with the number of elements and 
the Spearman correlation coefficient. The black dashed line cor- 
responds to the relation for resolved galaxies when assuming a 
solar metallicity (Eq. 10). The absence of metallicity estimates 
for a number of galaxies prevents us from using a metallicity- 
dependent Xqo factor. The red dots represent the 4 galaxies in 
sample A, with the bars indicating the full dynamical range 
probed by each of them. 



limited size of the sample, for most of the galaxies the relation 
derived for resolved galaxies provides a reasonable estimate of 
the face-on optical depth. The relations derived for both sam- 
ples prove to be very similar. This suggests that both resolved 
and unresolved galaxies follow relatively similar laws. Had they 
been markedly different, the derived relations would have shown 
obvious differences. This is especially important as sample B is 
made of only 27 galaxies and the scatter is large. A larger sample 
would be particularly useful to ascertain this result. Ultimately, 
the similarity between the relations suggests that the relations 
derived from resolved galaxies can be reasonably applied to un- 
resolved galaxies. 

To compare the dynamical range covered by resolved galax- 
ies versus unresolved ones, we have computed the global face- 
on optical depth for each galaxy in sample A. To do so, 1) we 
correct each pixel for the FUV attenuation using the estimate 
from CIGALE, 2) we sum the attenuation corrected flux in each 
pixel, 3) we compute the global FUV attenuation of each galaxy 
comparing the summed corrected and uncorrected FUV fluxes, 
and finally 4) assuming a slab geometry we compute the face-on 
optical depth. The red dots in Fig. 10 represent these integrated 
values and the bars around them the dynamical range covered by 
individual pixels within each of these galaxies. We see that our 
sample is well within the envelope described by sample B. Note 
that the mean gas surface densities for the 4 galaxies in sample 
A determined this way is much more secure than the values de- 
rived for sample B as they do not depend on assumptions on the 
distribution of the gas in galaxies. 
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7. Practical methods to compute the attenuation 

In this section we present some methods to compute the attenu- 
ation at all wavelengths in galaxies from the relations presented 
above, at kpc scales for resolved galaxies and for unresolved 
galaxies as a whole. The computation of the attenuation can be 
broken into several individual steps which we describe in detail 
hereafter: 

1. The mean gas surface density must be computed for each 
element (either an unresolved galaxy or an individual subre- 
gion) in the sample, taking into account both the atomic and 
molecular gas. 

(a) If either component is missing, it can be compensated for 
by assuming a certain molecular fraction. 

(b) In the case of unresolved galaxies, the computation of 
the gas surface density will be especially dependent on 
the actual distribution of the gas in the galaxy because 
this determines the area it should be averaged upon. 
We recommend averaging over 3 scale-lengths as we 
mentioned in Sect. 3.2.2: rm — 3 x 0.61 x r^s (Bigiel 
& Blitz 2012), and r Hl _ = 3x 0.2 x r 25 (Leroy et al. 
2008; Lisenfeld et al. 2011), with r 2 5 the optical radius 
at 25 mag arcsec~ 2 . 

2. From the gas surface density, the attenuation in the FUV 
band should be computed using one of the relations provided 
in this paper. 

(a) If an estimate of the metallicity is available, we strongly 
recommend the use of a relation explicitly dependent on 
the metallicity (Eq. 13). This is important as the metallic- 
ity affects the gas-to-dust mass ratio which has a direct 
consequence on the optical depth. We would like to stress 
as already mentioned in Sect. 5.3 that for consistency, 
the oxygen abundance estimator of Zaritsky et al. (1994) 
has been adopted throughout this paper. It is known for 
yielding rather high values compared to other estimators. 
In case estimates of the metallicity have been carried out 
with another estimator, we recommend either to recom- 
pute it with the Zaritsky et al. (1994) estimator, or to con- 
vert it with the formulas published in Kewley & Ellison 
(2008). 

(b) If no estimate of the metallicity is available, we recom- 
mend the use of a relation parameterised only on the gas 
surface density (Eq. 12). We caution that these relations 
have been derived on local galaxies that are relatively 
metal rich on average, 8.8 < 12 + logO/H < 9.4 with 
the Zaritsky et al. (1994) estimator, which could lead to 
large biases in the estimate of the optical depth if applied 
without correction to lower metallicity objects. 

3. From the FUV face-on optical depth, the FUV attenuation 
can be computed assuming a given geometry (Eq. 8 for a 
slab geometry). 

(a) As the relations we have derived between the gas surface 
density and the face-on optical depth depend on the ge- 
ometry we initially assumed, it is of utmost importance 
to adopt the same geometry 2 . 

(b) For consistency, we recommend the use the value of the 
albedo published by Draine (2003b). 

4. Finally, from the attenuation in the FUV band, it is possible 
to compute the attenuation at any wavelength assuming an 
attenuation curve. In case of a starburst galaxy, the starburst 
attenuation curve (Calzetti et al. 1994, 2000) can be applied. 

2 If another geometry is adopted, please contact the lead author of 
this paper to obtain updated relations applicable to this geometry. 



For high redshift galaxies, the recently derived attenuation 
curve of Buat et al. (201 la) can also be used. 

8. Comparison with relations from the literature 

In this section we examine how the relations we have obtained 
compare to some of the relations published in the literature. 

In the Milky Way it is possible to directly measure the quan- 
tity of gas towards individual objects and the extinction under- 
gone by these objects. In this case, relations between the opti- 
cal depth and the gas surface density take the form of a linear 
relation between these 2 quantities, and have an explicit depen- 
dency on the attenuation curve and the gas-to-dust mass ratio. 
Guiderdoni & Rocca-Volmerange (1987) expressed this as: 

1 A A A v E(B - V) (ZY 

Tj = X X X X X Ih, (16) 

1.086 A v E(B-V) N H \Z g J h v ; 

with ~E H the face-on gas surface density, E(B - V)/Nh the quan- 
tification of the reddening, which is generally calibrated in the 
Milky Way with a screen geometry, Z the metallicity, Z the so- 
lar metallicity, and s a wavelength-dependent parameter to take 
into account the variation of the shape and the amplitude of the 
reddening with the metallicity. Following Guiderdoni & Rocca- 
Volmerange (1987), we take s = 1.35 in the FUV band. If we 
assume an extinction curve with Ry = 3.1 and E(B - V)/Nh = 
1.81 x 10~ 2 mag (M pc -2 ) -1 from Giiver & Ozel (2009), in the 

V band and for a solar metallicity it reduces to: 

r y = 0.052 xT. H , (17) 

which is close to the widely used relation of Guiderdoni & 
Rocca-Volmerange (1987) and Devriendt et al. (1999): 

Ty — 0.059 x (18) 

We have to note that at very high column density, Ry may be 
larger as can be seen in dense molecular clouds. For instance 
in the Perseus cloud, Foster et al. (2013) measured an increase 
from R v - 3 for Ay - 2, to Ry - 5 for Ay = 10. At such 
attenuations, no significant quantity of UV photon could escape. 
However, on the scale of the current study, it should concern only 
a small fraction of UV photons. Otherwise it should be noted 
that CIGALE fits would fail catastrophically, which is not the 
case. We therefore think that the adopted value of R v = 3.1 is 
reasonable under the present conditions. 

Using such relations directly has several drawbacks to trace 
star formation in galaxies. First of all they are normalised to the 

V band which is a poor tracer of star formation. Computing the 
optical depth in the FUV, which is critical to compute attenu- 
ated and non-attenuated SFR, requires not only to assume an 
extinction curve (typically Milky Way, Large Magellanic Cloud, 
or Small Magellanic Cloud depending on the metallicity of the 
object), but also to take into account the mixing of different pop- 
ulations which are differently affected by the attenuation. Also, 
these relations are calibrated on individual lines of sight in the 
Milky Way or in the Magellanic Clouds assuming a screen ge- 
ometry, but not on large scale regions or even unresolved galax- 
ies that have more complex stars-dust geometries. 

To compare these relations to the ones we have obtained at 
kpc scales, we assume an extinction curve with tfuv/tv = 2.5 
which approximately corresponds to a Milky Way-like extinc- 
tion curve. Steeper extinction curves such as the LMC or the 
SMC would yield larger FUV optical depths. We present the re- 
sults in Fig. 1 1 . There are several striking differences that must 
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Fig. 11. Comparison of the relations between the face-on op- 
tical depth and the gas surface density obtained in this paper 
(solid lines, Eq. 1 2 and 1 3), with the relation from Guiderdoni 
& Rocca-Volmerange (1987) (dotted line), and with the relation 
from Giiver & Ozel (2009) (dash-dot line). The coloured lines 
show metallicity-dependent relations according to the oxygen 
abundance scale show on the right. For clarity, we have plot- 
ted the Giiver & Ozel (2009) relation only in the case of a solar 
metallicity. 

be noted. First of all, as discussed earlier at low gas surface 
density, the relations derived in this paper are asymptotically 
constant in the typical range 1.5 < Tfuv 5= 4.5. The relations 
from Guiderdoni & Rocca-Volmerange (1987) or Giiver & Ozel 
(2009) on the contrary fall linearly to t> yy — as 2# — > 0. This 
comes from the fact that these relations are calibrated along a set 
of individual line-of-sights in the Local Group, with an actual 
geometry that is close to the assumed screen geometry, with- 
out any mixing between different regions. Conversely the rela- 
tions derived in this paper are based on large regions in galaxies 
and naturally, undergo a strong mixing between high and low 
gas surface density regions that have a different face-on optical 
depth (Sect. 5.2). The difference at low column density shows 
in particular the expected effect of mixing compared to classi- 
cally used relations. In any case, at zero metallicity Tfuv = 
for all relations. The second feature we can note is that as the 
gas surface density increases, the metallicity-dependent rela- 
tion of Guiderdoni & Rocca-Volmerange (1987) yields, for the 
same metallicity, higher optical depths. Once again, the empiri- 
cal relations obtained in this paper more naturally take into ac- 
count the complexity of the galaxies, which we see makes dust 
less efficient at attenuating stellar radiation than the relation of 
Guiderdoni & Rocca-Volmerange (1987). 

9. Impact on galaxy formation and evolution 

The study of the LFs is one of the main tools to understand 
galaxy formation and evolution, either through observations or 
models. It is therefore important to understand the exact impact 
of the new relations we have derived in this paper on LFs output 
from SAMs. As mentioned in Sect. 1, UV LF obtained from hy- 
drodynamical simulations can be sensitive to the way the attenu- 
ation is computed (Finlator et al. 2006; Fontanot et al. 2009b) To 
understand the impact of different gas-attenuation prescriptions 
on models, we examine here how such relations would impact 



UV and IR LFs derived from cosmological simulations, as these 
are among the main targets to apply such relations. To do so we 
exploit the catalogue of galaxies of Guo et al. (2011) which is 
based on the Millennium (Springel et al. 2005) and Millennium- 
II (Boylan-Kolchin et al. 2009) simulations. 

9.1. Sample and observational parameters 

Before computing LFs we need 1) to select a sample of galax- 
ies from the simulated catalogue, and 2) derive their observa- 
tional parameters such as the mean gas surface density or the 
UV and IR luminosities. We have to stress that the objective 
is to test the relative effect of different attenuation-gas surface 
density relations but not to constitute highly accurate LFs in an 
absolute way, therefore the UV and IR luminosities will be built 
from relations widely used in the literature, that we will describe 
hereafter. The reproducibility of local LFs is discussed in ap- 
pendix H. 

As a first step, we select all star-forming galaxies (SFR > 0) 
at z = from the Guo et al. (2011) catalogue, which yields a 
sample of 708005 galaxies out of a full catalogue of 13191859 
galaxies. 

The Guo et al. (201 1) catalogue provides for each galaxy the 
cold gas mass in units of 10 10 M Q /h, and the size of the gas disk 
in units of Mpc//z, which is defined as 3 times the scale length, 
assuming an exponential disk, which is similar to the way we 
have determined the disk sizes for sample B. We assume h - 
0.73 to follow the parameters used for the simulation (Boylan- 
Kolchin et al. 2009). We compute the mean gas surface density 
simply by dividing the cold gas mass by the disk area. 

To compute the attenuation in the FUV band, we use the rela- 
tions between the face-on optical depth and the gas surface den- 
sity determined for sample A. The inclination of each galaxy is 
drawn randomly from a uniform distribution between 0° and 90°. 
As the catalogue of Guo et al. (2011) does not directly provide 
an estimate of the metallicity, we compute the oxygen abundance 
using the stellar mass-metallicity relation of Kewley & Ellison 
(2008) adapted for the Zaritsky et al. (1994) scale. 

Ideally, the FUV luminosity could be determined from the 
SFH of each galaxy. Unfortunately, the sampling in redshift is 
too sparse to trace accurately the variations of the SFH over 
the sensitivity timescale of the FUV. As a consequence, to com- 
pute the FUV luminosity, we assume that the FUV contamina- 
tion from previous star formation episodes is small and therefore 
that it traces the current star formation. Following Murphy et al. 
(2011), if we assume that galaxies have been forming stars at a 
constant rate over the last 100 Myr and follow a Kroupa (2001) 
IMF, we find: 

L(FUV) = 5.89 x 10 9 x SFR, (19) 

with L(FUV) the FUV luminosity in L G , and the SFR in 
M Q yr _1 . The use of another relation, such as the one of Salim 
et al. (2007) for instance, would yield slightly different FUV lu- 
minosities but would not affect the results. The IR luminosity 
can be determined from a relation between the UV-to-IR lumi- 
nosity ratio (IRX), and the attenuation. Several of such relations 
have been published based on local star-forming galaxies. They 
yield very similar results (Boquien et al. 2012). We adopt the 
formula determined by Buat et al. (2005): 

A FUV = -0.0333x 3 + 0.3522x 2 + 1.1960x + 0.4967, (20) 

with x = log L(IR)/L(FUV). To retrieve the value of L(IR), we 
simply invert this equation numerically. Galaxies being mostly 
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optically thin in the mid- and far-IR, the IR emission does not 
depend on the inclination of the galaxy. We therefore determine 
the value of L(IR) for an A F uv calculated for a face-on observer. 

9.2. Impact of relations between the face-on optical depth 
and the gas surface density 

To evaluate the impact of different relations between the face-on 
optical depth and the gas surface density, we compare in Fig. 12 
the FUV and IR LFs derived for the relation we derived for sam- 
ple A (Eq. 12 and 13) to the one presented in Guiderdoni & 
Rocca-Volmerange (1987); Devriendt et al. (1999). Note that in 
this section we select only galaxies that have 8.5 < logM* < 
11 because the stellar mass-metallicity relation of Kewley & 
Ellison (2008) is valid only in this range, reducing the catalogue 
to 54716 objects. This does not impact our results as we only 
examine here the relative variations on LFs induced by different 
relations. 

We note several striking features. First of all, the use of a 
metallicity-dependent estimator can make a large difference. For 
instance when selecting bright galaxies that have M(FUV) < 
-18, taking the metallicity into account systematically lowers 
the counts in each bin by -0.19 + 0.22 dex on average for the 
relation of Guiderdoni & Rocca-Volmerange (1987); Devriendt 
et al. (1999) and increases by 0.21 + 0.14 dex for our relations 
(Eq. 12 and 13). For fainter galaxies with M(FUV) > -18, 
even if the number counts are left nearly unchanged, respec- 
tively -0.00 ± 0.01 dex and -0.00 ± 0.08 dex, there is a change 
of behaviour around M(FUV) = -14.5. Fainter than this thresh- 
old the number counts are lower with the metallicity-dependent 
relation and higher when brighter than this threshold. A dif- 
ference in the number counts between metallicity-independent 
and metallicity-dependent relations is seen for the IR LFs. This 
shows the importance of taking the metallicity explicitly into ac- 
count, especially when studying the bright end of the LFs. When 
considering different metallicity-dependent relations between 
the gas surface density and the attenuation, we see that there 
are again significant differences that affect the overall LFs. The 
peak difference for UV bright galaxies reaches nearly 0.6 dex. 
Conversely in the IR, the difference in number counts remains 
small. 

To ensure that these results are not affected by the 0.2 mag 
width of the bins (smaller bins are more affected by random sam- 
pling effects), we have also examined the FUV LFs using a bin 
width of 0.5 mag. This size is similar to that of the observed FUV 
LF published by Wyder et al. (2005). We find that the previous 
results remain unaffected by an increase of the FUV bin width 
to 0.5 mag. 

These results are in line with those obtained by Fontanot 
et al. (2009b) and show the significant impact the relation be- 
tween the gas surface density and the face-on optical depth can 
have on the predicted LFs with peak differences in counts reach- 
ing about 1 dex. 

10. Summary and conclusions 

In this paper we have studied a sample of 4 resolved galaxies, 
and a sample of 27 unresolved galaxies in order to determine 
new but still simple relations between the gas surface density 
and the face-on optical depth at kpc scales, while also taking 
into account the effect of the metallicity. 

To do so, we have determined the gas mass surface density 
from HI and CO observations. In parallel, we have determined 
the face-on optical depth in two steps. 
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Fig. 12. Top: FUV LFs derived following the method outlined 
in Sect. 9.1. The blacks lines represent the LFs when using the 
gas surface density versus attenuation relation of Guiderdoni & 
Rocca-Volmerange (1987); Devriendt et al. (1999) whereas the 
blue lines represent the LFs when using the relations presented 
in the current paper, assuming a slab geometry. For both rela- 
tions, the solid line represent a relation parameterised only on 
the gas surface density whereas the dashed lines represent rela- 
tions paramaterised both on the gas mass surface density and the 
oxygen abundance. Bottom: IR LF, with the red lines represent- 
ing the use of the relations derived in this paper. 



1. We have determined the attenuation for each pixel in the 
resolved galaxies and for each unresolved galaxy using the 
CIGALE SED fitting code, which is based on an energy bal- 
ance principle between the attenuation in the UV-optical do- 
main and the emission of the dust in the FIR. 

2. From the attenuation, we have determined the face-on opti- 
cal depth assuming a slab geometry. 

We have then determined new relations between the gas sur- 
face density, and the face-on optical depth in resolved galax- 
ies, taking into account the metallicity indirectly by way of 
only metallicity-dependent Xco and gas-to-dust mass ratios, or 
also adding an explicit parameterisation. We have found that for 
galaxies with an inclination no larger than 50°, the assumptions 
made on the geometry have little impact on the derived attenua- 
tion. For larger inclinations, deviations over 30% can be found, 
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depending on the geometry. We can summarise our main find- 
ings as following: 

1 . The attenuation seems to primarily take place in regions that 
have a high molecular fraction, which is consistent with the 
fact that star-forming regions are more attenuated than more 
quiescent ones. 

2. The attenuation and the gas column density are affected dif- 
ferently by the averaging between different regions within 
a given resolution element in galaxies. This mixing effect 
strongly affects the relation between the face-on optical 
depth and the gas surface density, in particular in regions 
with a low surface density. This shows that oft-used rela- 
tions calibrated in the Milky Way are not adapted for use on 
large areas in galaxies or even entire galaxies. 

3. Unsurprisingly the impact of the metallicity is very signifi- 
cant, even for nearby spiral galaxies whose metallicity range 
is not extreme. As such it should be taken into account when- 
ever possible. 

4. Relations between the gas surface density and the face-on 
optical depth derived on resolved galaxies are able to repro- 
duce the trend for unresolved galaxies despite the scale dif- 
ference of these objects and the intrinsic difficulty to measure 
reliably the total gas surface density in unresolved galaxies. 

5. There is both a significant scatter around the best fits and 
galaxy-to-galaxy variations. This shows that there are ad- 
ditional physical parameters governing the relation between 
the face-on optical depth and the gas surface density. 

6. The new relations have a visible impact on LFs from 
the Millennium-II simulation when compared to the oft- 
used relations of Guiderdoni & Rocca-Volmerange (1987); 
Devriendt et al. (1999). This is especially visible on the 
bright end, and could therefore have a significant impact on 
the study of galaxy evolution through the use of large cos- 
mological simulations combined with SAMs. 

Finally, to compute the attenuation from the gas surface den- 
sity in practice, we recommend to apply the method laid out in 
Sect. 7. We remind that if another geometry is assumed, the lead 
authors of this paper should be contacted in order to obtain rela- 
tions adapted to this geometry. 
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Appendix A: The combination of FUV and dust 
emission as a proxy to the attenuation 

While we have shown in Sect. 3.1.2 that the determination of 
the FUV attenuation with CIGALE is reliable, we examine here 
whether and how the results we have obtained are affected by us- 
ing an alternative method to compute the attenuation and there- 
fore the face-on optical depth. Several relations exist in the liter- 
ature to quantify the attenuation from the FUV and the emission 
of the dust. These relations have generally been derived from 
the observation of unresolved star-forming galaxies (Cortese 
et al. 2008; Buat et al. 2011b; Hao et al. 2011). Boquien et al. 
(2012) have also derived such a relation but on subregions in 
star-forming galaxies using CIGALE. Even though it yields at- 
tenuations close to that of previous relations, we do not consider 
it here as it has also been determined using CIGALE. Among 
the relations from the literature, Boquien et al. (2012) found that 
the one provided by Hao et al. (2011) yields the largest differ- 
ence with the attenuation determined with CIGALE in resolved 
galaxies: (AAfuv) = 0.156 + 0.056 mag. We therefore adopt 
this relation as it will be indicative of the maximum deviation 
that could be induced by such an attenuation tracer compared 
to SED modelling. In Fig. A.l, we present the relations between 
the gas surface density and the FUV face-on optical depth deter- 
mined from the formula published in Hao et al. (2011), assum- 
ing a slab geometry. We see that using an alternative method to 
determine the attenuation leaves the relations qualitatively un- 
changed. After verification, the relations provided by Cortese 
et al. (2008) and Buat et al. (201 lb) show little qualitative differ- 
ence compared to when using the relation of Hao et al. (201 1). 
However the lower attenuations yielded by the relation of Hao 
et al. (2011) naturally affect the quantitative relations between 
the face-on optical depth and the gas surface density: 

T FUV = (2.353 ± 0.127) + (0.109 ± 0.009) x T, H , (A.l) 

and 

T FUV = [(1.859 + 0.095) + (0.029 ± 0.004) x 

x j q( 1 . 168±0.085)( [ 12+log O/H] -9.00) (A- 2) 

These relations yield smaller optical depths. Averaging over 
the entire sample, we obtain (Atfuv/tfuv) = 0.02 + 0.04 
when taking only the gas surface density into account, and 
(Atfuv/tfuv) - 0.10 + 0.07 and when taking into account both 
the gas surface density and the metallicity. The amplitude of 
these shifts is not larger than the error bars on the optical depths 
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Fig.A.l. Face-on optical depth in the FUV band determined 
from the FUV and dust emission (Hao et al. 2011) versus the 
atomic total gas mass surface density in units of M pc~ 2 . 



deduced from the uncertainties on the FUV attenuation provided 
by CIGALE. Therefore the use of different methods to determine 
the attenuation should only have a limited impact on the derived 
optical depth, with offsets of the order of 10%. 



Appendix B: Dust emission as a proxy to the gas 
surface density 

If the 21 cm emission is a direct tracer of the atomic gas, as 
mentioned earlier, the molecular gas is generally traced indi- 
rectly through CO emission. An alternative way to measure the 
total gas content in a galaxy is through the emission of the dust. 
Indeed, as gas and dust are intimately intertwined, assuming a 
given gas-to-dust mass ratio, one can retrieve the gas mass from 
the dust mass. Both methods have different biases. Whether the 
gas mass determination method affects our results or not is im- 
portant to estimate their reliability. 

One common method to determine the dust mass is to fit 
the FIR emission with a modified black body (e.g. Boselli et al. 
2002; James et al. 2002; Eales et al. 2010, 2012, and many oth- 
ers). It is described in Eq. B.l: 



fr 



1 



4ttD 2 



8tt/iv 3 
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with f v the flux density, D the distance of the object, h the 
Planck constant, v the frequency, c the speed of light in vac- 
uum, k the Boltzmann constant, T the temperature of the dust, k 
the dust emissivity taken at frequency v K , /3 the emissivity index, 
and Mdust the dust mass. Following Draine (2003a) we adopt 
k = 0.192 m 2 kg" 1 at v K = 857 GHz (350 fim). The value of 
jS is the subject of an intense debate in the literature. A value 
comprised between 1.5 and 2.0 is adequate for star-forming 
galaxies. Even though (3 can vary from region to region (Smith 
et al. 2012), Boselli et al. (2012) found that /3 = 1.5 best re- 
produces the FIR properties of the star-forming galaxies in the 
HRS sample. We therefore adopt their value. Note however that 
this emissivity has been determined from a dust model that as- 
sumes jS = 2 rather than /3 — 1.5, which may induce a bias in the 
dust mass. We determine the dust mass by fitting Eq. B.l to the 



data at 70 //m, 250 pm, and 350 /mi, using the curve_fit func- 
tion in the scipy library (Jones et al. 2001-), which is based on 
a Levenberg-Marquardt algorithm. This combination of bands 
allows us to obtain a good handle on both the dust temperature 
because these bands straddle the FIR emission peak, and on the 
dust mass as the 350 /mi band does not show a strong depen- 
dence upon temperature variations. In the case of NGC 2403, 
the 70 /mi emission comes from a distinctly different dust com- 
ponent rather than the dust emitting at longer wavelength. As a 
consequence, the dust mass within this galaxy will be underes- 
timated. In Fig. B.l, we compare the gas mass estimated from 
the dust mass, versus the one derived from HI and CO emission. 
In both cases we use the metallicity-dependent conversion fac- 
tors of Boselli et al. (2002) which assumes a gas-to-dust mass 
ratio of 160 for 12 + log O/H = 8.91. We see that there is a tight 
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Fig. B.l. Total gas surface density derived from a modified black 
body (Eq. B.l) fit of the FIR data from 70 /mi to 350 /mi versus 
the total gas surface density derived from the combination of 
HI and CO observations. The colour of each dot indicates the 
metallicity indicated by the scale on the right. The black lines 
show the 1-to-l relation between these quantities. 

correlation between the total gas surface density derived with a 
modified black body and from the combination of HI and CO 
observation as shown by the high value of the Spearman correla- 
tion coefficient for NGC 2403 (p = 0.79), NGC 4254 (p = 0.95), 
NGC 4321 (p = 0.96), and NGC 5194 (p = 0.95). We observe 
that the surface densities obtained using either method are sim- 
ilar even though there are some small offsets, generally towards 
lower masses except for NGC 2403. However they are small 
compared to the full dynamical range. These relations will be 
discussed in detail in Boquien et al. (2013, in prep.). 

To verify whether determining the total gas surface density 
affects the results, we plot in Fig. B.2 the relation between the 
face-on optical depth and the total gas surface density derived 
from a modified black body fit. When comparing to Fig. 7, we 
see that the results are very similar. The Spearman correlation 
coefficients for NGC 4254, NGC 4321, and NGC 5194 yield 
values close to what is found when combining CO and HI obser- 
vations. However, in the case of NGC 2403 there is a significant 
decrease from p = 0.72 to p = 0.45. This could be due to the 
lower metallicity of this galaxy ((12 + logO/H) = 8.90) com- 
pared to the rest of the sample ((12 + log O/H) - 9.20). As such 
there could be a significant quantity of "dark" molecular gas not 
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Fig. B.2. Face-on optical depth in the FUV band versus total gas 
surface density determined by fitting the FIR emission with a 
modified black body as described in Eq. B. 1. 



with P the angular power pattern of the beam, c/Q the solid angle, 
and mb denoting the main beam. For a gaussian beam with a full 
width half maximum size Aa and A6 along the major and minor 
axes, it easily comes that 



PdQ. = 

mb 4 In 2 



AaAd. 



Combining the 2 previous equations we get: 

SA 2 x21n2 

T B = 

knAaAS 

Inserting this relation into equation C.2: 
A 2 x 1.461 x 1(T 5 x21n2 



£h/ — 



knAaAS 



J>AV. 



(C4) 



(C5) 



(C6) 



Assuming that Aa and A6 are expressed in arcsec, 5, in Jy per 
beam, and AV in m s , a simple numerical application yields: 



8.85 
AaA5 



J>AV. 



(C7) 



traced by CO. Quantitatively, the relations we derive are close to 
Eq. 12 and 13: 

t fuv = (2.297 ± 0.059) + (0.092 ± 0.004) x I H , (B.2) 



and 



t fuv = [(2.206 ± 0.032) + (0.029 ± 0.002) x E H ] 

x j q(1.048±0.025)( [ 12+log O/H] -9.00) 



(B.3) 



These relations yield similar optical depths. Averaging over 
the entire sample, we obtain (Atfuv/ t fuv) - 0.09 + 0.05 
when taking only the gas surface density into account, and 
(Atfuv/tfuv) = -0.02 ± 0.06 and when taking into account 
both the gas surface density and the metallicity. The amplitude 
of these shifts is not larger than the error bars on the optical 
depths deduced from the uncertainties on the FUV attenuation 
provided by CIGALE. Therefore the use of different methods 
to determine the gas surface density should only have a limited 
impact on the derived face-on optical depth, with offsets of the 
order of 10%. 



Appendix C: Computation of the HI surface density 

In Sect. 3.2 we provide a relation to compute the HI mass sur- 
face density. It is derived in the following way. First we consider 
this well known relation linking Nhi, the HI column density in 
atoms ctrT 2 , to the brightness temperature Tg in K: 



Nhi = 1.823 X IO^Yt^AV, 



(CI) 



with AV in m s l . Considering that 1 M pc 2 = 1.248 x 
10 20 atoms cirT 2 , we have: 



E H / = 1-461 x 10~ 5 ^r B ,,AV. 



(C.2) 



The brightness temperature and the flux density per beam are 
related through the following equation: 



S = 



IkJn 



A- 



JOL- 

J J mb 



Pd£l, 



(C3) 



Appendix D: Impact of different CO transitions 

As explained earlier, to trace the molecular gas this study relies 
on the CO(2-l) transition. Yet different transitions may yield 
different results. Indeed the intensity of the emission of the var- 
ious CO transitions not only depends on the abundance of CO 
but also on the local physical conditions of the molecular gas, 
such as the density or the temperature. For a given quantity of 
CO, a high level transition will trace more closely the denser 
and warmer part of the gas. Disentangling between these effects 
is particularly challenging and would require a dedicated study 
of the sample with appropriate data. Here we therefore do not at- 
tempt to separate these parameters. We rather examine whether 
and how using different transitions would affect our results. 

To complement our CO(2-l) data, we need additional tran- 
sitions at a resolution that is similar or better than 30". This 
is a particularly limiting factor because CO(1-0) observations 
have a resolution twice as coarse as that of CO(2-l) for a given 
instrument. For higher order transitions, such as CO(3-2), ob- 
servations are difficult to carry out as they require optimal at- 
mospheric conditions. Nevertheless, 3 galaxies in our sample 
have adequate CO(3-2) maps to further our study. As part of 
the JCMT Nearby Galaxies Legacy Survey (Wilson et al. 2009), 
Bendo et al. (2010) published a CO(3-2) map of NGC 2403, 
and Wilson et al. (2009) maps of NGC 4254 and NGC 4321. For 
NGC 5194, Koda et al. (201 1) have published a high resolution 
single-dish CO(1-0) map which we adopt. Additional CO(1-0) 
observations of several galaxies in the sample have been pub- 
lished by Kuno et al. (2007). Unfortunately, their 15" beam is 
not Nyquist-sampled which makes the maps not suitable for a 
pixel-by-pixel analysis. In Fig. D. 1 we present the molecular 
gas surface density derived from the baseline CO(2-l) obser- 
vations and the CO(3-2) ones (CO(1-0) for NGC 5194). Note 
that as we divided our CO(2-l) intensities by a factor 0.89 to 
scale them to the CO(1-0) line. We also scale the CO(3-2) ob- 
servations to CO(1-0) by dividing the intensities by a factor 0.34 
(Wilson et al. 2009). There is visibly a strong correlation be- 
tween CO(2-l) and other CO lines, with the Spearman correla- 
tion coefficient ranging from p = 0.68 for NGC 2403 to p — 0.96 
for NGC 5 194. In Fig. D.2 we plot the relation between the face- 
on optical depth and the gas mass surface density. We see that 



20 



Boquien et al.: Towards understanding the relation between the gas and the attenuation in galaxies at kpc scales 





10 2 






a 




§ 




x-y) 


10 1 


o 




o 








W 






10- 


O 




a 




[M, 




x-y) 


10 1 


o 




u 




= 




w 





NGC 2403(CO 3-2} / 
N=894p = 0.68^/ 

Jw ■ 


NGC 4254(CO 3-21.//' 
N=432 p 0.9§gJ/ 


NGC 4321 (CO 3-2) . / 
N = 572 p = 0.894/ 


NGC 5194(CO l'-djuj^ - 
N=993 p=0.9jjp^ 



13: 



Quantitatively, the relations we derive are close to Eq. 12 and 



io u f 



10" 10 1 10 2 10" 10 1 10 2 

S tt (CO 2-1) [Mp pc~ 2 ] E a (CO 2-1) [M P pc~ 2 ] 

Fig. D.l. Molecular gas surface densities computed from CO(3- 
2) observations (CO(1-0) for NGC 5194) versus the molecular 
gas surface densities computed from the CO(2-l) transition. All 
CO lines have been scaled to that of the CO(1-0) in order to 
apply the same X co factor. The black line represent the 1-to-l 
relation where both lines yield the same molecular gas surface 
densities. 
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Fig. D.2. Same as for Fig. 3 with the molecular gas sur- 
face density determined from the CO(3-2) line (CO(1-0) for 
NGC 5194). 



qualitatively the use of different CO transitions does not affect 
the global trends between these quantities and that the qualita- 
tive relations remain unchanged. We note however that there are 
much fewer datapoints for NGC 5194, which makes the flatten- 
ing at low surface densities less evident. After inspection of the 
images, these are regions on the outskirts of the galaxy which 
are not in the field-of-view of the CO(1-0) map and have there- 
fore been eliminated from the analysis. We should note that as 
commented by Wilson et al. (2009); Bendo et al. (2010), there 
is a significant scatter in the ratio between the CO(3-2) line and 
lower CO transitions. It has been suggested that the CO(3-2) 
line are more closely associated with molecular gas fueling star 
formation than lower order lines, which may contribute to the 
pixel-to-pixel scatter. 



t fuv = (2.259 ± 0.063) + (0.099 ± 0.004) x (D.l) 



and 



tfuv = [(1-997 ± 0.033) + (0.033 + 0.002) x 2 H ] 

x j q( 1 .041 ±0.024)( [ 1 2+log O/H] -9.00) 



(D.2) 



These relations yield similar optical depths. Averaging over 
the entire sample, we obtain (Atfuv/tfuv) - 0.08 + 0.05 
when taking only the gas surface density into account, and 
(Atfuv/tfuv) = 0.05 + 0.04 and when taking into account 
both the gas surface density and the metallicity. The amplitude 
of these shifts is not larger than the error bars on the optical 
depths deduced from the uncertainties on the FUV attenuation 
provided by CIGALE. Therefore the use of different CO tran- 
sitions to determine the gas surface density should only have a 
limited impact on the derived face-on optical depth, with offsets 
of the order of 10%. 



Appendix E: Relation with the molecular fraction 

Whether the attenuation chiefly takes place in HI or H2 domi- 
nated regions gives us an interesting insight into radiation trans- 
fer in galaxies. In Fig. E.l we show the FUV attenuation versus 
the molecular fraction. The regions that have a low molecular 
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Fig. E.l. Relation between the molecular gas fraction and the at- 
tenuation. There is a clear trend, with increasing molecular frac- 
tion the attenuation increases, with the exception of NGC 2403 
which shows no trend. This is consistent with the trends ob- 
served in Fig. 3. The number of elements and the Spearman 
correlation coefficient are indicated at the bottom of each panel 
whereas the uncertainty on Apuv is shown in the bottom right 
comer. 

fraction generally correspond to peripheral regions. In the case 
of NGC 2403 we do not see any trend of the attenuation with 
the molecular fraction. Conversely the other 3 galaxies show 2 
regimes. Up to a molecular fraction of 0.6-0.8, there is little to 
no trend, the attenuation remaining constant. Beyond this limit, 
the attenuation strongly increases with the molecular fraction. 
This suggests that where the molecular fraction is high, the at- 
tenuation is strongly linked to the quantity of molecular gas. At 
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lower molecular fraction it is possible that the atomic gas plays 
a larger role in the attenuation. However, there could also be 
mixing effects if the FUV attenuation is driven by regions more 
closely linked to molecular gas even in low density regions. We 
will examine this possibility in more detail in Sect. 5.2. 

Appendix F: Impact of the geometry 

F.1. Attenuation-optical depth relations 

For this study, we have adopted a slab geometry. However, as 
mentioned in Sect. 5.1, other geometries can also be considered. 
First of all, we examine the case of a sandwich model. Following 
Boselli et al. (2003), the relation between the face-on optical 
depth t,\ and the attenuation is: 



A A = -2.5 log 



i-ft 



ft 



VI _ <^UTi/cos i 



(l + exp ( y/T— ojat a /cos iu 

(l - exp ( yi - co a t A /cos /)) 



(Fl) 



with f the height ratio between the dust and stars layers. Boselli 
et al. (2003) provide the following calibration based on observa- 
tions: 

ft = 1.0867 - 5.501 x 10 5 ^, (F.2) 

with A the wavelength in m. For the FUV band at 151.6 nm, f 
1.0033. When £ = 1, Eq. Fl is identical to the slab case (Eq. 8). 
As we limit the computation of the face-on optical depth to the 
FUV band (see Sect. 7 for the extension to other wavelengths), 
we will not discuss the sandwich geometry further. 

We also consider the case where dust is concentrated into 
clumps that are distributed following a Poisson distribution. In 
this case the relation between the attenuation and the optical 
depth is: 



A yl = -2.51og(exp[-n(l 



with n the average number of clouds along the line of sight, and 
t a c the optical depth of an individual cloud. 

In Fig. F. 1 , we compare the effective attenuation for various 
geometries, versus the mean optical depth t 1 ". At low mean opti- 
cal depths, all geometries yield the same attenuation. This means 
that the actual distribution of the dust, slab, clumpy, or even 
screen has little impact on the effective attenuation. However 
from an optical depth of unity, there is a strong divergence. For a 
slab geometry (solid black line in Fig. F. 1), A a cc log t"' whereas 
in the case of a clumpy distribution, A A °c h, when t" ! — > +oo. As 
a consequence in the latter case the resulting attenuation curve is 
more and more grey as the optical depth increases. As we see, in 
the case of a distribution of very opaque clumps, the maximum 
attenuation is directly determined by the number of clumps and 
not by the optical depth of individual clouds. 



F.2. Gas surface density and the optical depth for different 
geometries 

We compute the relation between the gas surface density and 
the face-on optical depth for screen, and clumpy geometries to 
understand the impact of the assumed geometries on these rela- 
tions. 

If we assume that all galaxies in the sample have the same 
average number of clumps along the line of sight, we need h > 3 
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Fig. F.l. Attenuation versus the mean optical depth t'" at a given 
wavelength A, for a face-on galaxy with a slab geometry (solid 
black line), and a clumpy geometry assuming a mean number 
of clouds on the line of sight of 1 (blue), 2 (green), 3 (red), or 
4 (cyan). A simple screen geometry is also shown for reference 
(dashed black line). For simplicity we have not taken into ac- 
count the effect of the albedo in this plot, that is we have assumed 
that co A = 0. 



to be able to reproduce the observations. This is slightly larger 
than what was found by Holwerda et al. (2007) for a sample of 
14 disk galaxies in the nearby Universe (n = 2.6). Considering 
h — 3 and h — 4, we find the following relations between the 
optical depth and the gas surface density: 



t fuv = (1.455 ± 0.022) + (0.065 ± 0.002) x Y. H , 



CF.3) and 



t fuv = [(1.361 ± 0.017) + (0.038 ± 0.001) x £ H ] 

x j q(0.448±0.0 14)( [12+log O/H]-9.00) 



(F.4) 



(F.5) 



for n = 3, and 

t fuv = (1.464 ± 0.018) + (0.052 ± 0.001) x !.„. (F.6) 

and 

t fuv = [(1.364 + 0.014) + (0.029 + 0.001) x 



x -y q(0.427±0.0 1 3)( [ 12+log O/H] -9.00) 



(F.7) 



for n = 4. 

When we compare to Eq. F.4 or Eq. F.6, Eq. 12 yields 
very different optical depths. For instance, if we assume Y, H - 
10 M pc~ 2 , Tfuv — 3.426 for a slab geometry and tfuv = 1-988 
for a clumpy geometry with h = 4, reaching Tfuv - 15.715 and 
t fuv = 6.704 respectively for I. H = 100 M pc~ 2 . This exem- 
plifies strong the model-dependence of such relations. In order 
to avoid biases, it is therefore particularly important to use the 
same geometry when computing the attenuation as the one that 
was assumed for the relation between the gas surface density and 
the optical depth. 
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Appendix G: Inversion of the relation between the 
attenuation and the face-on optical depth for a 
slab geometry 

In Sect. 5. 1 we introduced the Lambert W function to express the 
face-on optical depth as a function of the attenuation (Eq. 9). 
We indicate here how we have derived this equation and how 
the Lambert W function has been introduced. First of all, let's 
consider the equation we want to invert: 



A a = -2.5 log 



' 1 - exp (— VI - <^>at a /cos z'V 



VI - co at a /cos i 
Let r' = V 1 _ to at a /cos i. We easily get: 
= 1 - l(r a4A V, 



10 0M * x e 



-10' 
-10° 



,0.41, x e -T' + \0» M A x „-l(J 



T 

,0.4A, 



10 ( 



aAAi 



= r' - 10' 



,0.4.4, 



x e - 100 "' =(r'-10 0M ')x^' 



-10°- 4 4j 



(G.l) 



(G.2) 
(G.3) 
(G.4) 
(G.5) 



Using the W Lambert function, which gives the invert relation of 
f(w) = we w , this is equivalent to: 

r' - 10 a4A -' = w(-10 0AA i x e- 1004 " 1 ) . (G.6) 

Redeveloping t' and rearranging slightly we obtain the final re- 
lation: 



VI - <*>A 



[ W (-10 ^x e - 10 °^) + 10 " 4 ^]. (G.7) 



Appendix H: On the reproducibility of local LFs 

As stated before, the aim of this exercise is not to reproduce ex- 
actly the LFs of nearby galaxies but to examine the impact of dif- 
ferent gas surface density attenuation relations. Indeed, several 
intrinsic aspects of the simulations tend to incur biases affecting 
the LFs. 

A first well known problem is that semi-analytic relations 
can yield too much gas compared to what is observed in the lo- 
cal universe. To correct for this effect, following Obreschkow 
et al. (2009), we rescale the gas mass in each galaxy so that the 
total cold gas mass density at z — is consistent with the value 
derived by Obreschkow & Rawlings (2009): Q.f a s s = 6.0 x 10~ 4 , 
with Q"** normalised to the critical density for closure: 



Peril 



3H 2 
8^G' 



(H.l) 



H being the Hubble constant, and G the gravitational constant. 
For H = 73 km s" 1 Mpc _1 ,p cri( = 1.0 X 10~ 26 kg trT 3 . We then 
rescale the cold gas masses from the simulation by multiplying 
them with the factor f~ , with £ defined as: 



P 



gas 



(H.2) 



with p™"" the mean gas mass density of the simulation over the 
entire volume. We find that p™"' = 9.9 x 10" 30 kg trT 3 , lead- 
ing to £ = 1.5. This is very similar to the factor 1.45 found by 
Obreschkow et al. (2009) for the De Lucia & Blaizot (2007) 
catalogue which was based on the Millennium simulation. In 
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Fig. H.l. Mass function of the cold gas mass of the Millennium- 
II simulation (black lines), uncorrected (dashed line), and cor- 
rected (solid) for the mass excess. For comparison the observed 
cold gas mass function derived by Obreschkow & Rawlings 
(2009) is shown in red. All functions have been rescaled to 
h = 0.73, following the Millennium-II simulation parameters 
(Boylan-Kolchin et al. 2009). 



Fig. H.l, we show the Millennium-II cold gas mass luminos- 
ity. At high mass, the model slightly underestimates the cold 
gas mass function compared to the observed on by Obreschkow 
& Rawlings (2009). Conversely there is a very strong excess 
of galaxies with a lower mass of gas, up to about 1 dex for 
log M(H) = 7.0. 

The SFR distribution function, which is not affected by the 
attenuation, can also be compared to the observed ones in the 
local Universe. In Fig. H.2, we show the observed SFR distri- 
bution functions of Martin et al. (2005), Buat et al. (2007), and 
Bothwell et al. (2011). It turns out that the simulated SFR dis- 
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Fig. H.2. SFR function of the Millennium-II simulation (black 
line). For comparison, the Martin et al. (2005) function (cyan), 
the IR (red) and UV (blue) selected Buat et al. (2007) functions, 
and the Bothwell et al. (2011) functions observed in the local 
universe have also been plotted. All functions have been scaled 
to h = 0.73 and to a Chabrier (2003) IMF. 
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tribution function does not convincingly reproduce any of them. 
At the same time, they also have significant differences from an- 
other. The difference is the smallest at relatively high SFR, but 
the difference rapidly increases under SFR=10 M yr _1 , at the 
level of the "knee" of the distribution. For a SFR < 1 Mo yr _1 , 
the simulated SFR distribution function presents a much flatter 
slope than observed ones. At the same time, Guo et al. (201 1) re- 
marked that the Millennium-II simulation produces a large frac- 
tion of non-star forming dwarf galaxies compared to what is 
observed. This dearth of star-forming galaxies combined with 
the effect of the SFH - dwarf galaxies tend to have a more 
bursty SFH than large spirals which form stars at a steady rate 
- could explain such a flattening of the slope. At optical and 
near-infrared wavelengths which are increasingly less sensitive 
to the exact SFH, this problem has a much weaker impact. Guo 
et al. (201 1) showed that the Millennium-II simulations can re- 
produce convincingly the SDSS LFs derived by Blanton et al. 
(2005) in the g', r', i', and z' bands. However, in the star- 
formation tracing bands, these departures from the observed 
SFR distribution functions indicate that because we derive the 
extinction-free UV luminosity of each galaxy from the SFR as 
we saw in Sect. 9.1, we cannot compare directly our simulated 
UV and IR LFs with the observed ones. This is the reason why 
we do not attempt to carry out a comparison between our LFs 
and the observed ones. 
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